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Sluimary 

The  stabilization  of  boundary  layers  via  active  control  promises 
drag  reduction  with  relatively  small  power  requirements  compeured  to 
passive  boundary  layer  control  techniques.  To  date,  all  experimental 
active  control  investigations  have  concentrated  on  two-dimensional 
disturbances  in  the  laminar  boundary  layer.  To  further  improve  the 
level  of  disturbance  attenuation  in  the  boundary  layer  requires  am 
inherently  three-dimensional  approach.  The  goal  of  this  study  was  to 
control  the  boundary  layer  response  to  random  three-dimensional 
disturbances  introduced  near  the  leading  edge  of  a  flat  plate.  This 
goal  was  not  fully  achieved  in  the  shortened  time  frame  of  this 
investigation . 

The  basic  flow  structure  to  be  controlled  is  the  three- 
dimensional  wave  packet,  i.e.  the  boundary  layer  response  to 
localized  pulse  excitation.  Since  the  control  is  to  take  place  in  the 
linear  region  of  the  transition  zone,  any  conceivable  flow 
disturbance  can  be  synthesized  and  hence  cancelled  by  appropriate 
wave  packet  superposition. 

The  project  was  divided  into  two  main  p2urts,  (i)  numerical 
modelling  of  the  control  process,  and  (ii>  experimental  investigation 
of  the  most  promising  control  arrangement  identified  in  (i).  The 
numerical  modelling  was  an  integral  part  of  our  approach  and  allowed 
us  to  a  priori  assess  the  best  possible  control  performance 
achievable  with  a  given  design.  The  numerical  model  allows  rapid 
testing  of  different  spanwise  and  streamwise  distributed  detector  and 
actuator  arrays.  Using  this  model,  a  promising  control  configuration 
was  identified,  implemented  in  ^e  Cambridge  University  low 
turbulence  wind  tunnel  facility  and  tested  in  part  (ii)  of  this  work. 

The  experimental  investigation  achieved  signific2uit  control  of 
single  mode  disturbances.  In  fact,  control  to  about  twice  the  natural 
background  level  was  possible  for  frequencies  within  the  unsted>le  T-S 
band  (attenuation  up  to  32dB).  The  disturbemce  detection  was  based  on 
a  difference  microphone  arrangement  in  order  to  reject  acoustic 
disturbances  in  the  facility.  This  microphone  arrangement  proved  to 
be  ideally  suited  for  the  current  purpose.  The  actuators,  embedded 
miniature  loudspeakers,  proved  capable  of  controlling  moderate  sized 
insted^ility  waves. 

Wave  packet  control  was  attempted  but  failed  due  to  inadequate 
real  time  behaviour  of  the  electronic  filtering  circuit.  This  was 
traced  to  exceedingly  large  rise  and  settling  times  of  the  electrcxiic 
circuit . 
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1  .  Introduction 

The  possibility  of  "active"  control  of  instability  waves  in 
boundary  layers  has  attracted  the  renewed  attention  of  scientists 
since  the  early  80s.  Active  control  denotes  a  concept  in  which  the 
manifestations  of  flow  instability,  i.e.  the  Tollmien-Schlichting 
(T-S)  waves  are  influenced  by  direct  means  using  wave  superposition 
principles.  It  is  inherently  different  from  "passive"  methods  of 
control  which  modify  meaun  flow  features  to  achieve  steibility  with 
regard  to  small  disturbances. 

The  potential  benefit  of  active  rather  than  passive  control  of 
external  fluid  flows  (i.e.  boundary  layers)  lies  in  the  reduced 
power  requirement  of  the  former.  Even  though  passive  means  of  flow 
control  have  been  successful  in  flight  experiments  in  delaying 
transition  by,  say,  boundaury  layer  suction  or  wall  cooling  (air), 
the  undeniable  benefit  in  drag  reduction  for  practical  applications 
is  partly  negated  by  the  power  requirement  for  the  specific 
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technique  employed  (Antonatos,  1966;  VAiites,  Sudderth,  euid  Wheldon, 
1966;  Pfenninger  and  Reed,  1966;  Nenni  and  Gluyas,  1966;  Wagner  amd 
Fischer,  1983). 

The  power  requirement  for  an  active  control  configuration  is 
proportional  to  the  energy  content  of  the  flow  disturbances  amd 
thus  is  small  compared  to  a  passive  control  implementation. 
Nossenchuck  (1982)  increased  the  transition  length  in  a  flat  plate 
boundary  layer  in  water  by  25%  using  am  active  controller  with  1 
Watt  power  input.  The  same  delay  using  passive  control  (wall 
heating)  required  1900  Watts! 

Though  not  as  widely  known  as  the  passive  flow  control  methods 
mentioned  above,  active  flow  control  has  been  explored  since  the 
mid  60s  (Wehrmann,  1965;  Schilz,  1965/66)  and  re-emerged  in  the 
early  80s  (Milling,  1981;  Liepmamn,  Brown  and  Nosenchuck,  1982; 
Thomas,  1983;  Strykowski  amd  Sreenivasan,  1985;  Maestrello,  1985; 
Ladd  amd  Hendricks,  1988;  Ladd,  1990).  Improved  data-acquisition 
techniques  amd  advamces  in  micro-electrcmics  are  at  least  paurtially 
responsible  for  the  renewed  focus  on  these  techniques. 

The  early  experiments  by  Wehrmann  (1965)  and  Schilz  (1965/66) 
demonstrated  the  experimental  viability  of  active  control  by 
cancelling  two-dimensional  instability  waves  via  a  downstream  out- 
of -phase  controller  input  (through  flexible  wall  or  sound 
excitation).  Milling  (1981)  used  the  vibrating  ribbon  technique  to 
excite  T-S  waves  in  his  water  channel  experiment.  A  second  wire 
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downstream  introduced  an  "anti-wave”,  i.e.  a  wave  with  phase  and 
amplitude  adjusted  to  camcel  the  downstream  propagating  T-S  wave. 
Liepmann,  Brown  and  Nosenchuk  (1982)  used  a  novel  technique  to 
excite  T-S  waves  in  their  water  ttinnel  experiment.  Here,  flush 
mounted  surface  heaters  were  activated  to  excite  the  T-S  waves  in  a 
zero  pressure  gradient  flat  plate  boiandaury  layer.  Downstream 
cancellation  was  attempted  via  a  second  set  of  heater  strips  driven 
with  the  appropriate  anti-wave.  A  reduction  in  disturbance 
amplitude  and  a  delay  of  transition  was  achieved.  Liepmann  and 
Nosenchuck  (1982)  extended  the  technique  to  the  control  of  natxiral 
transition,  i.e.  the  T-S  waves  were  not  artificially  introduced  as 
in  the  above  mentioned  experimental  studies  but  were  generated 
naturally  by  the  background  free-stream  disturbcmce.  A  hot-film 
probe  was  used  as  a  sensor  downstreeun  of  the  heater  strip  which 
served  as  the  control  actuator.  The  sensor  output  was  analyzed  and 
a  signal  was  synthesized  to  drive  the  actuator.  This  feedback 
control  system  showed  T-S  wave  attenuation  or  reinforcement 
depending  on  the  relative  phase  of  sensor  20x3  actuator  signal. 
Thomas  (1983)  studied  the  influence  of  vibrating  ribbon  control 
(two-dimensional)  on  the  evolution  of  the  three-dimensional 
structures  in  the  transition  region.  He  showed  that  even  though 
transition  could  be  delayed  with  two-dimensional  active  control, 
the  weeik  three-dimensional  background  disturbance  field  interacted 
with  the  remainder  of  the  two-dimensional  primary  field  to  restart 
the  tremsition  process.  A  modified  vibrating  ribbon  technique  was 
used  in  a  wind  tunnel  experiment  by  Strykowski  and  Sreenivasan 
(1985)  to  excite  and  then  cancel  the  T-S  wave.  In  this  experiment. 
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the  ribbons  where  supported  in  slots  in  the  sxirface  of  the  plate 
and  hence  created  disturbances  more  like  localized  suction  and 
blowing.  The  experiment  was  successful  in  controlling  the  TS-wave. 
They  also  tried  to  use  the  heater  strip  technique  of  Liepmeuin  et  al . 
to  generate  and  cancel  T-S  waves  in  air.  However,  they  were  unable 
to  observe  any  flow  perturbations  using  this  technique.  This 
suggests  that  the  heater  strip  technique  might  not  be  suitaible  for 
active  control  actuation  in  zero  pressure  gradient  airflow 
applications.  In  boundary  layers  with  streamwise  pressure 
gradients,  however,  this  might  be  different.  In  fact,  Maestrello 
(1985)  has  used  heating  strips  to  generate  instability  waves  in 
boundary  layers  with  favorable  (1)  pressure  gradients.  He  also 
used  an  external  sound  field  to  generate  an  out-of-phase  control 
signal  which  led  to  T-S  wave  reduction.  The  recent  experimental 
investigations  by  Ladd  and  Hendricks  (1988)  and  Ladd  (1990)  are 
concerned  with  the  active  control  of  artificially  excited  and 
natural  occuring  instability  waves  on  em  axissymetric  body.  The 
experiments  were  carried  out  in  a  water  tunnel  using  hot-film 
detector  and  exciter  arrays.  A  feedback  loop  using  an  adaptive 
control  algorithm  was  used  to  successfully  control  two-dimensional 
disturbamces .  The  authors  discuss  the  strong  three-dimensional 
(azimuthal)  disturbance  components  and  envision  improved  results 
if  a  three-dimensional  control  scheme  were  implemented. 

Numerical  simulations  of  active  boundary  layer  control  have  been 
performed  by  various  investigators.  Most  of  them  integrate  the 
Navier-Stokes  equations  numerically  in  one  form  or  another.  A 
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review  [see  Krai  (1988)]  of  these  is  outside  the  scope  of  this  work, 
it  suffices  to  mention  that  until  recently  the  majority  dealt  only 
with  two-dimensional  disturbances  and  essentially  modelled 
vibrating  ribbon  experiments  [but  see  three-dimensional  numerical 
simulations  of  active  control  in  boundary  layers  by  Laurien  euid 
Kleiser  (1985),  Laurien  (1985)  and  Zang  and  Hussaini  (1985a, 
1985b)].  All  studies  confirm  the  feasibility  of  active  control 
transition  delay  or  enhancement. 


2.  Scope  of  Investigation 

The  above  survey  highligths  the  achievements  auid  limitations  of 
current  active  control  boundcUTr  layer  research.  One  serious 
limitation  of  all  the  reported  experimental  investigations  is  their 
restriction  to  two-dimensional  active  control.  This  limitation  is 
a  heritage  of  Schubauer  and  Skreunstad's  (1947)  experimental 
verification  of  the  nearly  two-dimensional  nature  of  T-S  waves  in 
the  early  stage  of  the  transition  process.  By  now,  however,  it  is 
well  established  that  the  transition  process  is  inherently  three- 
dimensional  and  the  concentration  on  two-dimensional  waves  seems 
>jn justified.  In  fact  the  above  cited  studies  by  Thomas  (1983)  and 
Ladd  (1990)  discuss  the  need  of  controlling  the  three-dimensional 
structures  in  order  to  advance  the  active  control  concept.  The 
classical  vibrating  ribbon  cannot  generate  a  three-dimensional 
control  input  and  hence  a  different  configuration  is  required.  The 
technique  of  flush  mounted  surface  heaters  used  by  Liepmann  et  al . 


-  6  - 


is  in  principle  suited  to  three-dimensional  control  actuation, 
however,  the  relative  large  size  of  the  surface  heaters  and  the 
rather  limited  power  input  possible  might  preclude  their 
successful  use  in  air  flow  applications. 

The  current  investigation  is  based  on  analytical  and  experimental 
work  on  three-dimensional  wave  packets  in  the  laminar  boundary 
layer  by  Caster  (1968,  1975,  1981,  1985).  Starting  conceptually 
from  a  local  three-dimensional  disturbcmce  in  the  boundary  layer 
produces  an  approach  more  suited  to  the  control  problem  at  hand. 
Since  it  is  possible  (in  linear  theory)  to  superpose  point  sources 
to  synthesize  arbitrary  forcing  fields,  the  pulse  response  can  be 
viewed  as  a  generic  flow  disturbance.  A  localized  pulse  generates 
a  wave  packet  which  propagates  downstream  in  the  boundary  layer 
amd  may  form  into  a  turbulent  spot.  Isolated  turbulent  spots  grow 
and  merge  to  form  a  fully  developed  turbulent  boundary  layer. 
Caster  (1975)  has  shown  that  the  wave  packet  development  in  its 
early  stage  can  be  modelled  accurately  by  linear  mode 
superposition.  Hence,  extending  the  ideas  of  T-S  wave  cancellation 
via  an  anti-wave,  the  cancellation  of  the  evolving  wave  packet  by 
superposition  with  an  "2mti-packet"  should  be  possible.  Viewing 
the  linear  wave  packet  as  fundamental  building  block  of  the 
transition  process  we  propose  the  use  of  wave  packet  superposition 
to  delay  boundary  layer  transition  due  to  aurtificially  excited  and 
-  ultimately  -  randomly  excited  disturbances . 

In  order  to  establish  the  feasibility  of  three-dimensional  random 
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disturbance  control,  the  following  issues  were  investigated  in  the 
course  of  this  project:  (a)  sensor/actuator  selection  for  control 
problem,  (b)  sensor /actuator  arrauigement  based  on  numerical  model 
of  control  system  including  boundary  layer  response,  (c) 
disturbance  detection  in  a  "noisy"  environment.  Finally,  a  wind 
tunnel  experiment  was  carried  out  in  order  to  evaluate  the  selected 
control  arrzmgement. 

This  report  is  divided  into  three  sections.  In  Section  3  below  the 
computer  model  of  the  linear  three-dimensional  control  process  it 
described.  Once  a  promising  control  configuration  was  identified 
through  the  modelling  process,  the  required  hardware  was  built  and 
an  experimental  investigation  was  launched.  In  Section  4  this 
experimental  investigation  is  described  in  detail.  Final,  in 
Section  5  conclusions  are  drawn  from  the  experience  gained  during 
this  investigation  and  a  course  of  future  reseeirch  is  suggested. 


3 .  Computer  Model  of  Active  Control  Process 

in  order  to  minimize  the  amount  of  experimental  "trial  and  error" 
necessary  to  find  a  viable  arrangement  of  local  surface  detectors 
emd  control  actuators,  full  use  was  made  of  computer  simulation  to 
test  ideas  and  optimize  (up  to  a  certain  point)  the  setup  before  the 
active  control  experiment  was  performed.  In  Subsection  3.1,  the 
specific  flow  control  configuration  will  be  presented.  The  control 
tas)c  will  be  discussed.  Sxibsectic^  3.2  deals  with  the  model  of  the 
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flow  field  response  based  on  the  three-dimensional  Orr-Sommerfeld 
stability  equation.  In  Subsection  3.3  we  digress  and  compare  our 
numerical  flow  field  response  with  the  experimental  flow  field 
response  obtained  in  a  preliminary  wind  tunnel  investigation.  In 
Subsection  3.4  we  use  the  model  of  our  flow  field  response  to 
tackle  the  boundaury  layer  control  problem  by  modelling  various 
spanwise  detector  amd  actuator  arramgements .  This  modelling  leads 
to  the  definition  of  the  most  "suitad)le"  control  arrangement  for 
the  control  experiment  described  in  Section  4.  (Listings  of  the 
computer  programmes  developed  in  the  course  of  this  project  caui  be 
found  in  Appendix  E.) 


3.1  Control  Task  Definition 

At  the  outset  it  is  necessary  to  specify  the  system  to  be  modelled 
and  controlled.  As  stated  adx>ve,  the  unstable  T-S  disturbances 
propagating  in  a  flat  plate  Blasius  boundairy  layer  are  to  be 
controlled.  Figure  1  shows  a  schematic  of  the  flat  plate  boundaury 
layer  and  introduces  the  control  problem. 

In  typical  flow  situations,  disturbamces  which  aure  naturally 
present  in  the  flow  environment  will  enter  the  boundary  layer  and 
generate  T-S  instability  waves.  The  processes  by  which  these 
external  disturbances  enter  the  boundary  layer  -  the  "receptivity" 
problem  -  aure  still  not  fully  understood  amd  urgently  require 
further  research.  Our  current  understanding  suggests,  however. 
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that  -  in  the  absence  of  any  wall  roughness  or  wall  humps  -  the 
external  disturbances  enter  the  boundary  layer  predominantly  in 
the  leading  edge  region.  In  order  to  model  this  process  properly,  a 
distributed  disturbamce  generation  in  the  region  downstream  of  the 
leading  edge  would  be  necessary.  For  the  current  investigation, 
however,  the  disturbance  generation  process  is  assumed  to  take 
place  at  a  certain  fixed  distauice  downstream  of  the  leading  edge. 
It  is  important  to  emphazise  that  the  artificial  disturbance 
generation  downstream  of  the  leading  edge  is  an  aid  in  the 
modelling  amd  experimental  phase  of  the  current  investigation. 
Ultimately  the  flow  disturbances  should  be  generated  by  "natural" 
means,  i.e.  turbulence  grids  upstream  of  the  plate. 

Once  the  disturbances  have  entered  the  boundary  layer,  they  will  be 
damped  or  amplified  as  they  propagate  downstream  (see  Subsection 
3.2)  and  will  reach  the  detector  location  in  the  flow.  Here  a 
spanwise  array  of  detectors  will  measure  the  disturbances  present 
in  the  flow.  In  essence,  the  continuous  time-dependent  2md 
spanwise  varying  flow  field  will  be  "discretized"  with  a  finite 
number  of  spanwise  detectors.  (Since  an  analogue  control  circuit 
will  be  used  no  time  discretization  takes  place.)  After  detection 
of  the  signal  it  will  be  analyzed  and  filtered  appropriately.  A 
control  signal  will  be  synthesized  and  injected  into  the  flow  at 
discrete  spanwise  actuator  locations. 

The  controlled  flow  field  will  be  filtered  by  the  boundaxy  layer 
and  the  control  performance  will  be  monitored  downstream  of  the 
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control  location  in  the  "target"  zone.  It  is  conceivable  that  a 
separate  row  of  detectors  in  this  target  zone  provides  the 
necessary  information  for  true  closed  loop  feedback.  In  the 
present  study,  however,  signals  obtained  at  the  target  location  are 
only  used  to  monitor  the  control  performance  and  are  not  fed  back 
into  the  control  circuit. 

The  equivalent  of  the  just  described  physical  system  is  represented 
in  block  diagram  form  in  Fig.  2a.  Script  H  denotes  a  system 
trcinsfer  function  and  round  and  square  brackets  identify  continuous 
and  discrete  representations  respectively.  HI  (k,w)  and  H2(k,w) 
denote  the  boundary  layer  response  to  the  external  forcing 
upstream  and  downstreeun  of  the  control  location.  Hl2(k,6>)  is  the 
boundary  layer  response  in  between  the  physical  location  of  control 
detection  and  actuation.  The  bottom  branch  of  the  loop  in  the  block 
inci'uies  deteyuor  controller  emd  actuator  transfer 
functions.  Detector  and  actuator  transfer  functions  are  each  split 
into  two  di.stinct  parts,  the  first  (1)  taking  into  acco”nt  the 
discrete  spanwise  nature  of  the  detection  and  actuation,  the  second 
part  (2)  dealing  with  the  finite  frequency  response  of  the  physical 
components.  Included  in  the  figure  are  some  possible  control 
strategies,  namely  gain  scheduling  or  adaptive  control.  None  of 
those  however  will  be  considered  in  this  report. 

A  simplified  block  diagreun  is  presented  in  Fig.  2b.  Here,  the 
loop's  complexity  is  reduced  by  neglecting  the  disturbance 
modification  due  to  the  flow  field  inbetween  detector  emd  actuator. 
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This  simplification  is  justified  if  detector  and  actuator  are 
physically  close  together  as  will  be  the  case  in  the  configuration 
chosen  for  our  study.  Furthermore,  in  this  simplified  case,  the 
finite  frequency  response  of  the  detectors  amd  control  actuators  is 
neglected,  i.e.  these  components  are  presumend  to  have  am  infinite 
bandwidth  flat  frequency  response  characteristic.  This  assumption 
is  justified  for  the  current  investigation  since  embedded 
microphones  and  loudspeakers  were  choosen  for  detection  and 
actuation  respectively  and  these  respond  essentially  flat  in  the 
relevant  T-S  frequency  bauid  [150  Hz  to  400  Hz]. 


3.2  Boundary  Layer  Impulse  response 

The  modelling  of  the  control  circuit  described  in  the  previous 
subsection  requires  knowledge  of  the  boundary  layer  response,  i.e. 
the  system  transfer  functions  HI,  HI 2  and  H2  (Fig.  2).  Since  the 
early  stage  of  the  boundary  layer  transition  process  is  described 
by  lineau-  theory,  the  impulse  response  of  the  boundary  layer 
provides  the  required  information.  Linear  systems  theory  allows  to 
construct  the  flow  response  to  any  arbitr2ury  input  by  simple 
convolution  of  this  input  with  the  impulse  response  function. 

Since  eui  impulsive  input  is  modelled  as  a  Dirac  spike,  the  initial 
condition  for  the  numerical  model  of  the  boundary  layer  impulse 
response  consists  of  a  flat  spectrum  of  unit  amplitude.  The 
downstream  development  of  each  Individual  Fourier  component  caui  be 
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determined  from  classical  linear  stadsility  theory.  The  pulse 
response  at  a  certain  downstream  location  is  obtained  by 
superposing  the  contribution  of  all  Fourier  components  at  this 
point.  The  lineeu:  theory  describing  the  three-dimensional 
disturbance  development  in  a  flat  plate  boundary  layer  is  based  on 
the  Orr-Sommerfeld  equation  invoking  Squire's  transformation  to 
reduce  the  three-dimensional  problem  to  am  equivalent  two- 
dimensional  one  (Mack,  1969).  Solving  the  resulting  eigenvalue 
problem  provides  the  complex  wavenumber  "a"  as  a  function  of  local 
Reynoldsnumber  amd  forcing  frequency  tJ.  The  impulse  response  at  a 
certain  downstream  location  "xj^  as  a  function  of  spanwise  position 
is  thus  described  by 


u'  (time, 2-span) 


Here  u'  denotes  the  disturbance  velocity  component  parallel  to  the 
wall  in  mean  flow  direction.  Approximately  10000  eigenvalues  "a" 
need  to  be  calculated  in  order  to  represent  the  temporal  and 
spamwise  development  properly.  The  integral  in  the  exponent  is 
evaluated  using  a  Romberg  integration  procedure.  In  order  to 
rapidly  calculate  the  laorge  number  of  required  eigenvalues,  am 
efficient  algorithm  (Caster,  1978)  was  used  which  expresses  the 
eigenvalue  dispersion  relationship  in  a  complex  double  series, 

■  2  Z  -  l^lf  (ci'  *  b'  - 

n  m 
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llie  spatial  eigenvalue  evaluation  required  Newton-Raphson 
iteration.  Using  the  Caster  double  series,  the  10000  eigenvalues 
were  cdstained  in  2Lbout  15  minutes  on  a  SUN-Sparc  station.  It  should 
be  noted  that  only  eigenvalues  falling  within  a  certain  area 
surrounding  (2md  of  course  inside)  the  neutral  st2d>ility  curve  were 
considered  since  the  other  highly  damped  modes  do  not  contribute 
significantly  to  the  impulse  response  far  away  from  the  source. 

Using  the  described  procedure,  the  required  bound2iry  layer  impulse 
response  functions  HI  and  H2  were  calculated.  Figure  3  shows 
contour  plots  of  the  spanwise  distribution  of  the  disturbance  u- 
velocity.  Also  shown  is  the  convolution  H1*H2,  corresponding  to 
the  flow  response  at  the  target  location  (see  Fig.  1 )  due  to  pulse 
forcing  at  the  exciter  location.  Hie  modelled  boundaury  layer 
impulse  response  produces  the  expected  wave  packet  shape  (Caster 
and  Crant,  1975;  Caster;  1975).  In  the  following  subsection,  this 
predicted  wave  packet  shape  will  be  compared  more  closely  with  the 
one  obtained  from  experiment. 


3.3  Comparison  of  Model  and  Experiment 

The  above  described  boundary  layer  response  calculation  forms  the 
building  block  for  the  modelling  of  the  control  arrangement. 
Hence,  before  starting  extensiv  simulations  based  on  this  model  it 
is  sensitive  to  compare  the  predicted  impulse  response  of  the 
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boundary  layer  flow  with  suitable  experinents .  To  this  end  a  test 
plate  was  built  and  experiments  were  performed  in  the  Cambridge 
University  low  turbulence  reseeurch  wind  tunnel.  Details  of  the 
hardware  are  described  in  Appendix  C  and  the  experiment  itself  is 
presented  in  Appendix  D.  The  experiment  was  set  up  to  model  the 
parameters  of  our  numerical  model,  i.e.  disturbance  generaticxi  at 
Re^f  =  870,  and  measuring  position  at  1235,  corresponding 

roughly  to  the  fictitious  "control"  location  of  Re^^  1260. 
Measxirements  of  the  fluctuating  u-velocity  component  where 
obtained  with  a  constant  temperature  hot-wire  positioned  just 
outside  the  boundary  layer  where  the  T-S  eigenfunctions  exhibit  an 
outer  (flat)  maxima.  In  the  experiment,  the  flow  was  excited  by  a 
short  duration  (0.5ms)  computer  generated  pulse  via  an  embedded 
loudspeaker.  The  results  are  presented  as  isometric  views  and 
contour  plots  for  the  u-velocity  fluctuation. 

Figure  4  shows  the  experimental  wave  packet  results  from  ensemble 
averaged  records.  The  equivalent  flow  field  calciilation  based  on 
our  numerical  model  is  shown  in  Fig.  5.  Since  the  model  based  on 
the  Orr-Sommerfeld  eigenvalue  prc^lem  does  not  provide  absolute 
amplitude  levels,  the  amplitude  level  in  the  simulation  was 
adjusted  to  the  experimental  maximum.  Comparing  experimental  and 
numerical  result  shows  remarkable  agreement.  In  particulaur  the 
location  of  the  wave  packet  within  the  time  record  is  very  well 
reproduced  by  the  model.  ITiis  can  be  seen  even  more  clearly  in  Fig. 
6  which  shows  a  comparison  of  experimental  result  and  numerical 
result  on  the  packet  centre  line.  Considering  that  the  model  is 
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based  on  a  locally  parallel  boundeiry  layer  approximation,  the 
agreement  is  excellent. 


3.4  Control  Simulation 

Having  established  that  the  numerical  model  of  the  flow  response 
and  the  actual  flow  response  obtained  by  wind  tunnel  experiment  are 
nearly  identical,  the  simulation  of  the  control  arrangement 
introduced  in  Subsection  3.1  can  be  performed. 


3.4.1  General  Performance  Considerations 

Conceptually,  the  control  circuit  (Fig.  1 )  consists  of  the  following 
components  (in  flow  direction):  (i)  event  detectors,  (ii)  ccMitrol 
actuators,  emd  (iii)  residue  disturbance  detectors  which  are 
located  in  the  "target"  zone  downstream  of  the  actuators.  The 
event  detectors  upstream  of  the  actuators  are  required  since  the 
disturbance  wave  packets  (events)  appear  at  random  tines.  We  aire 

* 

searching  for  viable  control  geometries,  i.e.,  we  use  the  numerical 
model  to  rapidly  explore  the  advantages  and  disadvantages  of 
different  streamwise  and  spanwise  detector /actuator  arrai^ements . 
Ttie  numerical  "control"  process  proceeds  as  follows:  the  flow 
response  at  the  detector  location  (Re^»*  1235)  due  to  any  tine 
dependent  input  disturbance  f(t)  is  obtained  by  convoluting  f(t) 
with  the  apropriate  system  transfer  function  (HI ) .  As  canonical 
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case,  we  choose  the  pulse  as  input  disturbance,  i.e.  f(t)  «  S (t).  At 
the  control  locatiOT  we  "detect"  the  disturbance  and  can  "control" 
by  arbitrarily  modifying  the  numerical  signal  at  different  spanwise 
locations.  This  "controlled"  disturbance  is  then  convoluted  with 
the  sytem  transfer  function  H2  in  order  to  provide  the  disturbance 
flow  in  the  target  zone  at  R^.*  1780  (Fig.  1 ).  m  this  first  attempt 
to  esteUblish  a  viable  control  geometry  the  individual  detector  and 
actuator  dynamics  is  not  t20cen  into  account.  Hence,  "ideal"  control 
is  simulated,  i.e.  the  components  of  the  control  circuit  are 
assumed  to  have  an  undistorted,  infinitely  wide  frequency  response 
with  zero  response  delay  (see  Fig.  2b).  Three  different  open- loop 
control  configurations  are  explored: 

-  Case  1,  the  ideal  controller  is  located  at  the  centre-line, 

-  Case  2,  in  addition  to  the  controller  at  the  centre-line,  two 
controllers  are  placed  symmetrically  12mm  to  each  side  of  the 
centre-line,  and 

-  Case  3,  a  total  of  four  controllers  are  used  which  are  placed 
symmetrically  at  6mm  and  21mffl  above  and  below  the  centre-line,  m 
order  to  assess  the  control  performance,  the  pseudo-energy  (u'  )  in 
the  target  zone  is  integrated  over  the  whole  disturbance  flow 
field.  The  cost  function  in  the  target  zone  is  minimized  by 
"optimal"  choice  of  controller  gain  and  phase.  Without  any  control. 
Pig.  7a  shows  the  modelled  flow  respcxise  at  the  target  location. 
iTte  pseudo  disturbance  energy  in  the  field  is  15.6.  This  serves  as 
the  refererice  case  to  assess  the  ccmtrol  performance.  Figure  7b 
shows  the  best  possible  control  for  case  1 ,  a  single  controller  on 
the  centre-line.  Even  though  the  fringes  of  the  packet  can  not  be 
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controlled,  the  disturbance  energy  is  reduced  by  a  factor  of  3  to 
5.0.  Figures  8a,  b  show  the  best  ccMitrol  possible  for  the  three 
controllers  and  four  controllers  configuration  respectively.  With 
four  controllers,  the  residual  disturbance  energy  is  reduced  to  0.6 
which  is  just  4%  of  the  uncontrolled  level.  The  sensitivity  of  our 
open  loop  control  to  v2iriation  in  amplifier  gain  and  controller 
phase  is  shown  in  Fig.  9.  As  expected,  proper  phase  control  is 
crucial  for  a  successful  active  control  implementaticxi.  It  is 
interesting  to  note  that  the  minimum  in  the  amplifier  gain  curve 
(Fig.  9a)  is  rather  shallow  Indicating  a  moderate  cost  penalty  for  a 
controller  gain  mismatch.  For  comparison,  the  dotted  line  in  Fig. 
9a,  b  indicates  the  control  sensitivity  for  the  case  of  a  purely 
two-dimensional  single  mode  wave  controlled  with  an  ideal  two- 
dimensional  control  arrangement.  Obviously,  a  control  amplifier 
gain  of  unity  allows  complete  cancellation  for  this  hypothetical 
case.  A  properly  matched  amplifier  gain  is  more  critical  for  this 
case  then  for  the  three-dimensional  wave  packet  cases  presented 
above,  while  the  phase  sensitivity  is  about  equal. 

From  the  simulations  described  in  this  subsection,  it  was  decided 
that  a  spanwise  spacing  of  about  10mm  would  provide  an  acceptable 
balance  of  achievable  control  performance  and  hardware  complexity. 
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3.4.2  Performance  of  chosen  Configuration 

The  "ideal"  ccmtrol  circuit  used  above  did  not  take  into  account  the 
physical  spacing  of  detectors  and  actuators.  It  was  based  on  the 
simplified  block  diagram  in  Fig.  2b.  Die  streamwise  spacing  of 
detector  and  actuator  array  can  be  included  into  our  numerical 
model.  Furthermore,  the  detecticxi  vas  modelled  as  an  ideal  point 
detection.  Bmbedded  microphones  were  choosen  to  provide 
disturbance  pressure  detection.  Experience  from  previous  own 
experiments  and  the  work  by  Kendall  (1990)  shows,  that  acoustic 
disturbances  present  in  the  wind  tunnel  environment  produce  wall 
pressure  fluctuations  an  order  of  magnitude  larger  than  those 
associatfid  with  the  instability  waves.  These  acoustic  pressure 
fluctuations  have  to  be  eliminated  from  our  signal  in  order  to 
allow  detection  of  the  relevant  disturbances,  the  instability 
waves .  This  was  done  using  two  micrc^ones  spaced  a  "small"  (on  an 
acoustic  scale)  distance  apart  in  streamwise  direction  and  switched 
such  that  the  signals  are  substracted  from  each  other.  The 
streamwise  distance  was  in  effect  chosen  to  be  roughly 
(0.5*lambda),  "lambda"  being  the  wavelength  of  the  most  amplified 
instability  wave  component  at  the  detector  location.  Hence  this 
double  microphone  zurangement  will  effectively  reject  the  Icxig 
acoustic  wavelength  but  will  enhance  the  detectability  of 
disturbances  with  wavelengths  in  the  relevant  instability  wave 
range.  The  behaviour  to  be  expected  for  this  arrangement  is  shown 
in  a  Bode  plot  in  Fig.  10.  The  frequency  for  which  the  difference 
microphcxie  arrangement  is  tuned  is  denoted  by  F,  and  a  6dB 
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amplification  is  achieved  for  this  particular  frequency.  For  lower 
frequencies  corresponding  to  longer  wavelengths,  the  eunplitude 
response  is  lower,  asymptotically  approaching  a  typical  first  order 
attenuation  of  20dB  per  frequency  decade.  It  is  conceptually 
helpfull  to  view  the  double  microphone  arrangement  from  a  finite 
difference  standpoint.  Essentially  the  difference  microphone 
approximates  the  derivative  of  the  passing  signal  at  the  midpoint 
between  the  micrqphone  location.  Integrating  the  difference 
microphone  signal  electronically  should  recover  the  exact  signal 
at  the  midpoint  (without  the  constant  contribution  from  the 
acoustic  contamination)  providing  the  "ideal"  response  indicated  in 
the  graph.  The  actual  combined  response  of  difference  microphones 
and  integrator  is  constructed  by  simple  addition  of  the  appropriate 
curves  in  the  Bode  plot.  It  should  be  noted  that  in  principle  the 
signal  detected  with  the  difference  microphone  arrangement  is 
distorted  by  folding  (aliasing)  due  to  higher  disturbance 
frequencies  with  shorter  wavelength  them  can  be  resolved  with  the 
selected  streamwise  separation  of  the  microphones.  However,  since 
the  unstable  frequency  band  is  rather  narrow  and  frequencies  which 
exceed  by  a  factor  of  2  or  higher  are  strongly  damped,  this 
spatial  aliasing  is  Insignificant. 

b  order  to  minimize  interference  between  actuator  loudspea)cer  and 
detector  microphone,  it  was  decided  to  place  the  actuator 
loudspeakers  staggered  in  the  spanwise  direction  with  respect  to 
the  micrc^hone  array  (Fig.  10).  It  is  then  sensible  to  synthesize 
the  control  actuator  input  by  a  weigthed  average  of  the  nearest 
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difference  microphone  detectors. 

Including  the  difference  microphone  with  integrator  and  the 
spanwise  staggering  of  the  detector/actuator  arrangement  in  our 
numerical  model  and  taking  furthermore  the  boundary  layer  flow 
development  in  the  controller  region  (impulse  response  function 
H12  in  Fig.  2a)  into  account,  the  performance  of  the  full  control 
system  was  analyzed.  Figure  11a  shows  the  uncontrolled  model  flow 
field  at  the  target  location  due  to  excitation  at  the  exciter 
location  near  the  leading  edge.  The  input  at  the  exciter  location 
was  chosen  as  a  random  function  in  time  and  spanwise  distribution 
in  order  to  present  a  more  realistic  flow  picture.  (Of  course, 
testing  the  control  performance  with  a  single  pulsed  point  source 
exciter  reveals  exactly  the  same  information  since  the  system  is 
modelled  linearly.)  As  before.  Fig.  11a  shows  a  contour  plot  of  the 
disturbance  u-velocity  distribution  as  a  function  of  time  emd 
spfuiwise  positicxi.  For  the  sake  of  clarity,  the  center  time  trace 
and  the  spanwise  distribution  at  a  particular  point  in  time  are 
extracted  and  plotted  separately  below  and  to  the  right  of  the 
contour  plot  respectively.  Considering  the  temporal  structure 
first,  the  signal  looks  plausible  in  comparison  to  typical 
oscilloscope  traces  of  natural  (i.e.  random)  transition 
experiments.  Concentrating  on  the  spanwise  distribution  it  is 
appaurent  that  the  structure  is  predominantly  two-dimensional  in 
nature,  i.e.  the  wave  fronts  are  aligned  almost  perpendicular  to 
the  mean  flow  direction.  This  again  is  expected  from  linear 
stability  theory  since  two-dimensional  disturbances  are  in 
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principle  the  most  amplified  ones.  However,  the  structure  shows 
sufficient  spanwise  modulation  to  xiiin  any  control  attempt 
based  solely  on  a  disturbance  detection  at  one  spanwise 
location.  In  other  words,  this  plot  confirms  the  necessity  of  the 
distributed  control  approach  chosen  in  the  current  investigation  if 
cancellation  of  the  disturbauice  structure  over  a  significant 
surface  area  is  to  be  successful. 

Ttie  same  flowfield  at  the  target  location  with  the  control  active 
is  shown  if  Fig.  11b.  Obviously  no  total  cancellation  is  achievable 
with  the  choosen  detector/actuator  arrangement.  However,  the 
maximum  amplitude  in  the  whole  flow  field  was  reduced  to  18%  of  its 
uncontrolled  value  and  the  overall  pseudo  disturbemce  energy  was 
reduced  to  about  4%  of  the  uncontrolled  result.  The  achieved 
control  has  to  be  viewed  as  the  best  possible  control  in  an  ideal 
setting  and  it  is  to  be  expected  that  the  implementation  in  a  wind 
tunnel  experiment  will  degrade  the  performeuice .  Nevertheless,  the 
predicted  amplitude  reducticwi  to  about  one-fifth  of  its 
uncontrolled  value  justifies  the  implementation  of  the  control 
system  in  a  wind  tunnel  experiemnt. 
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4 .  Active  Control  Experiment 

An  insert  for  the  test  plate  was  designed  and  built  which 
incorporated  the  detector  and  actuator  arrangement  arrived  at  by 
numerical  modelling  as  described  zdxsve.  The  mechanical  and 
electronic  hardware  used  is  described  in  more  detail  in  Appendix  C. 
The  experimental  investigation  was  divided  into  three  different 
stages  described  below,  m  Subsection  4.1,  the  suitability  of  the 
embedded  miniature  loudspeakers  for  the  control  of  instability 
waves  was  tested.  Subsection  4.2  deals  with  the  important  issue  of 
disturbance  detection  and  noise  rejection.  Implementation  of  the 
full  control  arrangement  including  both  microphone  detectors  and 
loudspeaker  actuators  is  described  in  Subsection  4.3. 


4 . 1  External  Control 

As  a  first  step  to  a  successful  control  of  boundary  layer 
instability  waves,  it  is  important  to  ascertain  that  the  control 
actuators,  i.e.  the  miniature  loudspeakers  embedded  in  the  plate 
can  indeed  produce  the  disturbance  level  required  for  c2uicellation 
at  the  actuator  locaticxi.  In  order  to  do  so,  an  "external"  control 
experiment  wais  devised,  bypasssing  disturbance  detection  entirely. 
A  continuous  single  mode  disturbance  of  known  fixed  frequency  was 
introduced  at  the  exciter  location.  At  the  control  location,  the 
detectors  were  uncoupled  from  the  circuit.  The  control  signal 
driving  the  control  loudspeakers  was  generated  externally  with  a 
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frequency  generator  set  to  the  same  frequency  as  the  exciter 
frequency.  Finally  the  phase  and  amplitude  of  this  external 
"controller"  was  manually  adjusted  until  disturbance  cancellation 
was  achieved  as  monitored  by  a  hot-wire  on  the  centre-line  at  the 
target  location.  This  procedure  was  repeated  for  different 
frequencies  and  exciter  amplitudes  in  order  to  establish  the 
maximum  disturbance  amplitude  which  could  be  cancelled  using  the 
chosen  actuator  arrangement.  It  transpired  that  the  chosen 
actuators  could  control  disturbances  of  moderate  amplitude.  This 
moderate  amplitude  disturbances  were  generated  by  driving  the 
exciter  speaker  neair  the  leading  edge  with  roughly  40%  of  the 
maximal  voltage  allowable  for  linear  distiurbance  generation.  In 
other  words,  driving  the  exciter  speaker  2.5  times  as  hard  produces 
a  disturbance  2.5  times  as  big  (hence  the  system  response  linearly) 
but  the  control  actuators  were  not  able  to  cancel  it  at  the 
controller  location  downstream. 

Figure  12a  shows  a  hot-wire  trace  of  an  uncontrolled  disturbance 
generated  by  simultaneously  feeding  three  discrete  frequencies 
into  the  exciter  speaker.  Prior  to  this  superposition,  each  of 
these  three  frequencies  had  been  controlled  individually  by 
external  control  as  described  above.  If  the  system  is  linear, 
switching  on  the  externally  tuned  controllers  of  all  three 
frequencies  simultaneously  should  eliminate  the  signal.  Figure  12b 
shows  the  hot-wire  trace  of  the  controlled  result  and  comparison 
with  the  unforced  reference  case  (not  shown)  reveals  that  the 
residual  disturbance  in  the  controlled  case  is  twice  the  size  of 
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the  natural  background  noise.  Hence  from  this  experiment  we  could 
conclude: 

(a)  the  control  actuators  are  suitable  for  control  of  moderatly 
sized  instability  waves,  and 

(b)  tb**  ^ysten  behaves  linearly  as  expected,  i.e.  linear 
superposition  holds. 

As  a  next  step,  the  spemwise  distributicxi  of  the  control  'Vredge" 
generated  by  a  finite  number  of  spanwise  actuators  was  examined  in. 
an  external  control  setting.  In  Fig.  13,  the  spanwise  distribution 
of  the  hot-wire  rms  integrated  over  the  relevant  T-S  is  shown. 
The  particular  case  represents  superpose<'  forcing  and  control  of 
two  distinct  frequencies.  The  external  controller  was  driving  four 
control  actuators  simultaneaously  (indicated  by  black  triangles  on 
the  e±>cissa),  only  half  of  the  flow  field  is  shown  in  the  graph.  In 
the  uncontrolled  case  we  see  the  expected  roll-off  in  tue  spanwise 
direction.  In  the  controlled  configuration,  the  very  local  nature 
of  the  control  near  the  centre-line  is  appeurent.  It  is  interesting 
to  note  that  there  is  no  performance  penalty  at  the  fringes  of  the 
controlled  regicxi  due  to  possible  in-phase  amplification  of  the 
flow  disturbcince . 


4.2  Disturbance  Detection 

The  quality  of  the  sensor  signal  is  of  crucial  importance  for  good 
control  performance.  This  requires  proper  sensing  of  disturbance 
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amplitude  and  phase  as  well  as  adequate  noise  rejection.  The  later 
point  is  particuleu'ly  important  due  to  the  dominant  acoustic 
disturbance  field  present  in  the  wind  tunnel  environment.  As 
described  above,  a  double  microphone  arrangement  was  chosen  as 
sensor,  the  small  (acoustic  scale)  st.eamwise  spacing  of  the 
microphones  allowing  rejection  of  acoustic  disturbances.  In  order 
to  ascertain  the  "coherence"  of  the  relevant  flow  disturbance,  i.e. 
T-S  instability  wave,  and  detector  signal,  a  hot-wire  was  placed  at 
the  centre  of  the  double  microphone  in  streamwise  position  but 
displaced  sideways  in  spau  by  8mm.  This  spanwise  displacement  was 
chosen  in  order  to  avoid  any  influence  of  the  local  flow  field 
generated  by  the  hot-wire  support  structure  on  the  microphone 
signal.  A  pulse  disturbance  was  generated  at  the  exciter  locaticxi 
near  the  leading  edge  and  three  channels  were  recorded 
simultaneously  providing  (i)  the  hot-wire  signal,  (ii)  the  upstream 
microphone  signal  of  the  double  microphone  arrangement,  and  (iii) 
the  downstream  microphone  signal  of  the  same  arrangement. 

Figures  14a  and  14b  show  a  single  realization  of  the  hot-wire  trace 
and  the  ensemble  average  of  64  realizations  of  the  hot-wire  signal 
respectively.  Ihe  hot-wire  trace  shows  a  wave  packet  riding  on  a 
strong  background  disturb^ulce  field  of  comparable  magnitude.  The 
generated  wave  packet  is  clearly  discernible  after  ensemble 
averaging  the  signal.  However,  a  single  hot-wire  at  the  detector 
location  would  not  be  suitable  for  disturbance  detection  due  to 
insufficient  signal  to  noise  ratio.  Figures  15a  and  15b  show  a 
single  realization  of  a  single  microphone  and  of  the  dotible 
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Microphone  arrangement  respectively.  This  data  were  taken 
simultaneously  with  the  single  shot  hot-wire  trace  shown  in  Fig. 
14a.  No  wave  packet  can  be  detected  in  the  single  microphone 
signal,  rendering  a  single  microphone  totally  iinsuit2d>le  for  our 
purpose  of  instability  wave  detection.  However,  using  the 
Microphone  in  conjunction  with  the  second  micro(dK»ie  placed 
downstream  such  that  the  difference  signal  is  recorded,  a  spatial 
filter  of  excellent  quality  is  obtained.  In  the  single  realization 
double  microphone  trace  shown  in  Fig.  15b,  the  instability  wave 
packet  clearly  rises  above  the  background  noise  level.  The  signal 
to  noise  ratio  is  similar  to  the  one  obtained  by  ensemble  averaging 
the  hot-wire  with  64(1)  realizations. 


4.3  Open  Loop  Control  Performauice 

Given  the  control  actuator  performance  discussed  in  Subsection  4.1 
and  the  signal  quality  of  a  single  shot  difference  microphone 
arrangement  discribed  in  the  previous  subsection,  the  active 
control  of  instability  waves  in  the  boundary  layer  should  be 
achievable  with  the  present  hardwau-e  configuration.  In  this 
subsection,  we  describe  the  control  performance  achieved  using  the 
described  arrangement.  We  concentrate  on  the  control  performaunce 
at  the  centre-line.  Due  to  time  cxxistraints,  no  spanwise  traverses 
were  recorded.  The  whole  system  thus  looks  as  follows  (see  Fig.  1 ). 
Excitation  takes  place  with  a  single  "point  source"  exciter  on  the 
plate  centre-line.  At  the  control  location  detection  via  one  set  of 
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embedded  double  microphones  on  the  centre-line  t2d(es  place.  iTtis 
signal  is  then  filtered  (see  Appendix  C  for  hardware  details)  and 
fed  to  four  control  actuators  simultzmeously  which  are 
symmetrically  distributed  across  the  centre-line  2uid  are  10mm 
apart.  The  control  performance  is  assessed  in  the  target  zone 
downstream  of  the  control  location  with  a  hot-wire  placed  above  the 
plate  centre-line  at  the  outer  edge  of  the  boundary  layer. 

In  the  first  set  of  experiments  single  mode  control  was  attempted, 
using  the  full  control  arrangement.  Figure  16a  and  b  show  the  hot¬ 
wire  trace  at  the  target  location  with  control  off  and  on 
respectively.  The  almost  complete  C2uicellation  of  the  disturbance 
signal  is  apparent.  In  fact,  taking  the  Fourier  transform  of  the 
controlled  signal  and  comparing  the  integrated  rms-value  over  the 
T-S  frequency  band  of  the  controlled  and  the  uncontrolled  case 
reveals  a  reduction  of  about  a  factor  of  10  in  disturbance 
amplitude,  i.e.,  a  factor  of  100  in  disturbance  energy.  Comparing 
the  controlled  case  with  a  reference  case  without  any  disturbance 
excitation  (not  shown)  reveals  that  the  control  achieves 
attenuation  of  the  disturbance  amplitude  to  about  twice  the  natural 
background  level.  It  is  thus  apparent  that  our  control  arrangement 
is  as  good  as  the  external  control  arrangement  described  in 
Subsection  4.1  in  eliminating  the  instability  wave.  It  should  be 
pointed  out,  that  in  order  to  achieve  this  good  control 
performance,  a  phase  shifter  within  the  ccxitrol  circuit  (see 
Appendix  C)  had  to  be  adjusted  quite  caurefully.  The  control  circuit 
was  tested  for  different  exciter  frequencies  auid  amplitudes  and 


-  28  - 


ccMitrol  to  twice  the  background  level  was  consistently  possible  if 
the  internal  phase  shifter  was  adjusted  to  a  slightly  different 
phase  lag  for  each  frequency.  However,  even  if  the  phase  shifter 
was  not  precisely  tuned,  i.e.,  it  was  set  to  a  const2uit  average 
value,  the  control  of  single  mode  frequencies  over  the  whole 
unstable  T-S  band  was  still  significant,  i.e.,  a  factor  of  4  in  u- 
disturbance  amplitude  for  each  frequency. 

In  carrying  out  these  experiments,  it  became  apparent  that  the 
proximity  of  disturbance  detector  2uid  actuator  leads  to 
undesirable  crosstalk,  i.e.,  the  two  control  actuators  nearest  the 
detector  location  influence  the  detector  reading.  The  extent  of 
interference  is  shown  in  Fig  .17.  Detector  output  as  a  function  of 
control  actuator  driving  voltage  is  shown.  In  the  absence  of  any 
cross-talk,  the  detector  output  should  be  independent  of  the 
actuator  driving  level.  In  Fig.  17a  the  interaction  is  appeirent. 
Using,  however,  two  control  actuators  further  apart  from  the 
detector  microphones.  Fig.  17b  shows  that  the  Interference  prc^lem 
is  eliminated.  Since  the  ccxitrol  performance  using  the  arrangement 
of  Fig.  17a  did  not  deteriorate  profoundly,  no  attempt  was  made  to 
eliminate  this  interference  problem  by  further  filtering.  Hie 
second  problem  which  was  unexpected  was  the  varying  phase  lag 
requirement  in  the  ccxitrol  circuit  for  different  frequency 
components.  From  the  beginning  it  was  expected  that  a  phase  shift 
would  be  required  to  take  account  of  the  local  receptivity 
behaviour  at  the  actuator  location.  Xn  other  words,  since  It  is  not 
known  what  the  (complex)  "transmission"  coefficient  is  between  a 
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pressure  di8turb2uice  generated  at  the  wall  and  the  downstream 
disturbance  u-veloclty  component  if  would  be  fortuitous  if  it  were 
identically  zero.  It  was  unexpected  that  this  transmission 
coefficient  is  a  function  of  frequency.  Careful  measurements  of 
this  transmission  coefficient  were  performed  such  that  a  proper 
f^ase  shifting  all-pass  filter  can  be  designed  at  a  later  stage. 

As  the  last  experiment  in  this  investigation  full  wave  packet 
control  was  attempted.  The  difference  between  controlling  single 
mode  disturbances  and  a  full  wave  packet  is  the  appe2a'2uice  of  the 
latter  at  a  certain  a  priori  unknown  point  in  time.  Hence,  the  time 
constant  of  the  electronic  ccxitrol  circuit  plays  a  significant 
role.  The  influence  of  this  time  constant  was  underestimated  in 
our  control  arrangement  and  led  to  the  failure  in  controlling  the 
wave  packet.  It  was  however  not  only  the  larger  than  expected  rise 
time  of  our  electronics  which  caused  this  failure.  In  euSdition,  the 
settling  time  was  exceeding  the  passing  time  of  the  wave  packet 
thus  causing  a  ringing  of  the  ccxitrol  actuator  which  intrcxiuced 
additicxial  disturbances  into  the  boundary  layer  flow .  The  severity 
of  this  problem  can  be  seen  from  Fig.  18.  The  top  trace  shows  the 
detector  signal  and  the  wave  packet  can  be  identified  in  the  trace. 
The  spikes  seen  in  this  trace  are  due  to  a  ground  cxxinection  which 
was  left  floating  by  mistake.  This  was  corrected  in  later  attempts 
and  should  be  ignored  for  the  current  argument.  The  bottom  trace 
shows  the  cxxitrol  actuator  signal  driving  the  embedded  cxmtrol 
speakers.  The  delayed  start  of  the  ctxitrol  signal  and  the  Icxig 
settling  time  after  the  signal  has  passed  are  apparent.  This  Icxig 
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settling  tine  was  traced  to  the  narrow  band  filtering  associated 
with  the  employed  phase  shifting  device  and  could  not  be  corrected 
in  the  tine  available  for  this  experiment.  Tlie  time  delay  due  to 
the  Bessel  filters  was  found  to  be  in  the  order  of  1msec  which  is 
too  large  for  the  system  to  be  able  to  control  frequencies  in  the 
neighbourhood  of  300Hz  (period  3msec8)  effectively.  It  is 
important  to  note,  that  the  failure  of  the  attempted  wave  packet 
control  is  not  based  on  any  unexpected  flow  behaviour  but  purely 
due  to  the  inadequate  electr^ic  circuit.  Fortunately,  having 
concluded  this  experiment  and  having  specified  the  circuit 
requirements  it  is  feasible  to  construct  an  improved  electronic 
control  circuit  taking  into  account  the  lessons  learned. 


5 .  Conclusions  and  Outlook 

The  stabilization  of  boundary  layers  via  active  control  promises 
drag  reduction  with  relatively  small  power  requirements  compared 
to  passive  boundary  layer  control  techniques.  To  date,  all 
experimental  active  control  investigations  have  concentrated  on 
two-dimensional  disturbances  in  the  laminar  boundary  layer. 
However,  to  further  improve  the  level  of  disturbance  attenuation  in 
the  boundary  layer  requires  an  inherently  three-dimensional 
approach.  Ihe  work  presented  in  this  report  addresses  this  issue, 
the  goal  of  this  study  was  to  control  the  boundary  layer  response 
to  random  three-dimenslcmal  disturbances  Introduced  near  the 
leading  edge  of  a  flat  plate,  this  goal  was  not  fully  achieved  in 
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the  shortened  tine  frame  of  this  lnvestigatl<^.  (The  original 
proposal  set  out  a  three  year  programme,  however,  two  years  were 
funded . ) 

The  basic  flow  structure  to  be  controlled  is  the  three-dimensional 
wave  packet,  i.e.  the  boundary  layer  response  to  localized  pulse 
excitation.  Since  the  control  is  to  take  place  in  the  linear  region 
of  the  transition  zone,  any  conceivable  flow  disturbance  can  be 
synthesized  and  hence  C2mcelled  by  appropriate  wave  packet 
superpositicn . 

The  project  was  divided  into  two  main  parts,  (i)  numerical  modelling 
of  the  control  process,  and  (ii)  experimental  investigation  of  the 
most  promising  control  arrangement  identified  in  (i).  The 
numerical  modelling  was  an  integral  part  of  our  approach  and 
allowed  us  to  a  priori  assess  the  best  possible  control  performance 
achievable  with  a  given  design.  The  numerical  model  allows  rapid 
testing  of  different  spanwise  and  streamwise  distributed  detector 
and  actuator  arrays.  Using  this  model,  a  promising  control 
ccffifiguration  was  identified,  implemented  in  the  Cambridge 
University  low  turbulence  wind  tunnel  facility  and  tested  in  part 
(ii)  of  this  work. 

The  experimental  investigation  achieved  significant  control  of 
single  mode  disturbances.  In  fact,  control  to  about  twice  the 
natural  background  level  was  possible  for  frequencies  within  the 
unstable  T-S  band.  The  disturbance  detection  was  based  on  a 
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difference  microphone  arrangement  in  order  to  reject  acoustic 
disturbances  in  the  facility.  This  micrc^>hone  arremgement  proved 
to  be  ideally  suited  for  the  current  purpose.  The  actuators, 
embedded  miniature  loudspeakers,  proved  capable  of  controling 
moderate  sized  instability  waves.  Slight  cross  talk  between 
actuator  and  detector  arrang^ent  was  observed  but  was  not  so 
severe  as  to  adversely  affect  the  control  perform2nce. 

Wave  packet  control  was  attempted  but  failed  due  to  inadequate  real 
time  behaviour  of  the  electronic  filtering  circuit  which  could  not 
be  rectified  in  the  experimental  time  remaining.  This  failure  was 
due  to  exceedingly  large  rise  and  settling  times  of  the  electronic 
circuit.  The  former  caused  the  ccmtrol  to  "kick  in"  after  about 
half  the  wave  packet  had  passed  the  actuator  location,  while  the 
latter  caused  excessive  "ringing"  of  the  actuator  after  the  wave 
ptacket  had  passed  thus  in  itself  creating  a  disturbance  in  the  flow 
field. 

It  should  be  emphasized  that  the  failure  to  control  the  wave  packet 
is  not  due  to  any  unforseen  flow  behaviour  but  solely  due  to 
inadequate  electronic  circuitry.  Redesigning  the  electronic 
filters  is  therefore  the  necessary  first  step  in  achieving  wave 
packet  control  locally  at  the  plate  centre-line.  Increasing  the 
physical  distance  between  control  detectors  and  actuators  is 
advisable  in  view  of  (a)  eliminating  cross  talk  and  (b)  allowing  for 
more  processing  time  in  the  control  loop.  After  control  at  the 
centre-line  is  achieved,  it  is  suggested  to  carefully  explore  the 
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spcinwise  distribution  amd  compare  this  with  the  numerical  model 
developed  in  the  course  of  this  project.  Afterwards,  several 
spanwise  detector/actuator  blocks  can  be  combined  to  provide 
complete  spanwise  control  of  distributed  random  disturbance 
fields. 
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Appendix  A:  Figures 
List  of  Figures 

Fig.  1  Schematic  of  boundary  layer  control  aurrauigement . 

Fig .  2  Block  diagram  representation  of  boundeoy  layer  control 
system . 

(a)  full  system 

(b)  simplified  basic  system 

Fig.  3  Simulaticxi  of  flow  response  due  to  point  source  excitation 
at  Re^>=  870.  Dashed  lines  correspond  to  negative  values. 
ST1  -  System  transfer  functicxi  from  Re^>-  670  to  1235. 
ST2  -  System  transfer  function  from  Re/o  1235  to  1780. 
ST1*ST2  -  Disturbance  flow  response  at  1780  due  to 

point  source  excitation  at  Rej<=  870  (STI 
convoluted  with  STII). 

Fig.  4  Streamwise  velocity  r' actuations  (u')  measured  at 
outer  edge  of  bov  ry  layer  (Rey=  1235). 

Single  Pulse  excitation  at  Rej<s  870. 

Fig.  5  Numerical  modelling  of  streamwise  velocity 

fluctuaticxi  (Re^'s  1235).  Single  Pulse  excitation 
at  Ref.  =  870 . 

Fig .  6  Comparison  of  disturbance  u-velocity  fluctuation  along 
wave  packet  centre-line. 

Fig.  7  Simulated  flow  response  (u')  in  target  zone  (Refi=  1780), 
dashed  lines  correspond  to  negative  values. 

(a)  no  control  (dash-boxed  insert  of  Fig.  3  rescaled) 

(b)  case  (1 ):  single  controller  on  centre-line,  contour 
levels  as  in  (a) 

Fig.  8  Simulated  flow  respcxise  (u'}  in  target  zone  (Ref>*  1780), 

dashed  lines  correspond  to  negative  values.  Contour  levels 
as  in  Fig.  10a.  Relative  size  of  wave  packet  at  control 
location  Re^**  1235  and  speuiwise  controller  positioning  is 
indicated  in  the  left  part  of  plot  area. 

(a)  case  (2):  three  controllers  configuration 

(b)  case  (3):  four  controllers  configuration. 

Fig.  9  Performance  sensitivity  to  controller  gain  and  phase 
mismatch  for  different  control  arrangements. 

(a)  gain  mismatch 

(b)  phase  mismatch 
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Fig.  10  Bode  plot  of  amplitude  and  phaise  behaviour  of  difference 
microphone  arrangement  includiikg  ideal  integrator. 

(a)  amplification  vs.  normalized  frequency 

(b)  ph^e  shift  vs.  normalized  frequency 

Fig.  11  Simulated  flow  response  (u')  in  target  zone  (Rejc=  1780), 
due  to  random  disturbance  excitation  in  spanwise  position 
and  time.  Dashed  lines  correspond  to  negative  values. 

Time  trace  along  centre-line  and  spanwise  distribution 
at  intermediate  time  extracted  to  bottom  and  right  of 
plot  respectively. 

(a)  uncontrolled  reference  case 

(b)  controlled  case 

Fig.  12  Centre-line  hot-wire  trace  at  target  location. 

Excitation  with  three  discrete  frequencies . 

(a)  uncontrolled  case 

(b)  "external"  control  on 

Fig .  1 3  Integrated  hot-wire  rms  in  T-S  frequency  band  at  target 
location  versus  spanwise  probe  position. 

"External"  control  of  two  superposed  frequencies. 
Actuator  location  at  control  location  indicated  by 
full  triangles  on  abcissa. 

Fig.  14  Hot-wire  trace  at  detector  location  (centred  in 

streamwise  position  between  double  microphone  but 
displaced  sideways  by  Smm  to  avoid  interference). 

Pulse  excitation  near  leading  edge. 

(a)  single  realization  hot-wire  record 

(b)  ensemble  average  of  64  realizations 

Fig.  15  Microi^one  trace  at  detector  location  centre-line. 
Pulse  excitation  near  leading  edge.  Simultaneously 
recorded  with  hot-wire  trace  in  Fig.  14a. 

(a)  Single  realization  of  single  micrc^phone 

(b)  Single  realizaticxi  of  difference  microphone 

Fig.  16  Centre-line  hot-wire  trace  at  target  location. 

Single  mode  excitation  near  leading  edge. 

(a)  uncontrolled  signal 

(b)  controlled  signal  using  control  loop. 
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Fig.  17  Actuator /detector  cross-talk.  Detector  output  as 
a  function  of  externally  driven  actuator  forcing 
level . 

(a)  two  actuators  symmetrical  to  detector  microphones, 
spanwise  displaced  by  Smm . 

(b)  two  actuators  symmetrical  to  detector  microphones, 
spanwise  displaced  by  15mm. 

Fig.  18  Settling  time  of  control  circuit. 

(a)  Detector  trace  showing  passing  of 

wave  packet.  [Spurious  spikes  in  signal  are 
due  to  floating  ground  connection  in  this 
particular  case  and  should  be  ignored.] 

(b)  Synthesized  control  actuator  driving  signal 
showing  long  settling  time. 

Fig.  Cl  Aluminium  test  plate  for  active  control  experiments. 

Fig.  C2  Cambridge  University  Low  Turbulence  Research  Tunnel. 

Fig.  C3  Electronic  hardware  for  active  control  experiment. 

Fig.  D1  Boundary  layer  response  in  packet  centre-line  at  Re  =  1380. 
Data  analogoue  filtered  between  20Hz  and  2000Hz. 

(a)  single  realizaticxi 

(b)  ensemble  average  of  189  realizations 

Fig.  D2  Power  spectral  density  versus  frequency  plot  of  ensemble 
averaged  data  shown  in  Fig.  Dib. 

Fig.  D3  Measured  u- fluctuations  at  outer  boundary  layer  edge  at 
RejjC!  1380.  Ensemble  averaged  (189)  and  digitally  filtered 
between  80Hz  and  500Hz . 

(a)  perspective  view 

(b)  contour  plot,  negative  contours  dashed 

Fig.  D4  Measured  u- fluctuations  in  the  boundary  layer  at  Re  >  1235 
for  twin  pulse  Interaction  experiment .  Two  pulses  of  same 
strength  buth  opposite  sign  excite  flow  at  Re^  880 . 
Ensemble  averaged  data  digitally  filtered  between  80Hz  and 
500Hz. 

(a)  perspective  view 

(b)  contour  plot,  negative  ccm tours  dashed 

Fig.  D5  Measured  u-fluctuations  in  the  boundary  layer  at  Ref=  1235 
for  triple  pulse  interaction  experiment.  Three  pulses  of 
same  strength  (3cm  apart  spanwise)  excite  flow  at  Re^  880. 
Middle  pulse  of  opposite  sign.  Ensemble  averaged  data 
digitally  filtered  between  80Hz  and  500Hz . 

(a)  perspective  view 

(b)  contour  plot,  negative  cot  tours  dashed 
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Fig.  1  Schematic  of  boundary  layer  control  arrangement. 


(a) 


Fig.  2  Block  diagram  representation  of  boundary  layer  control 
system. 

(a)  full  system 

(b)  simplified  basic  system 


Fig.  3  Simulation  of  flow  response  due  to  point  source  excitation 
at  Re^»=  870.  Dashed  lines  correspond  to  negative  values. 
ST1  -  System  transfer  functicxi  from  Re^«=  870  to  1235. 
ST2  -  System  transfer  functicm  from  Re/#s  1235  to  1780. 
Sti*st2  -  Disturbamce  flow  response  at  Re^>>  1780  due  to 
point  source  excitation  at  Re..«  870  (STI 
convoluted  with  STn).  • 


EXPERIMENTAL  RESULTS  -  SINGLE  SOURCE 


Fig.  5  Numerical  modelling  of  streamwise  velocity 

fluctuation  (Reio  1235).  Single  Pulse  excitation 


.  7  Simulated  flow  response  (u')  in  target  zone  1780) 

dashed  lines  correspond  to  negative  values.  * 

(a)  no  control  (dash-boxed  Insert  of  Fig.  3  rescaled) 

(b)  case  (1 ):  single  controller  on  centre-line,  contour 
levels  as  in  (a) 
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Fig.  8  Simulated  flow  response  (u')  in  target  zone  (R^»  1780), 

dashed  lines  correspond  to  negative  values.  Contour  levels 
as  in  Fig.  10a.  Relative  size  of  wave  paclcet  at  control 
location  Re^*«  1235  and  spanwise  controller  positioning  is 
indicated  in  the  left  part  of  plot  area. 

(a)  case  (2):  three  controllers  configuration 

(b)  case  (3):  four  controllers  configuration. 
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Fig.  10  Bode  plot  of  amplitude  and  {rfiase  behaviour  of  difference 
microphone  arrangement  including  ideal  integrator. 

(a)  amplificaticm  vs.  normalized  frequency 

(b)  phase  shift  vs.  normalized  frequency 


Fig,  11  Simulated  flow  response  (u')  In  target  zone  (Re^<«  1780), 
due  to  random  disturbance  excitation  In  spanwlse  poslticxi 
and  time.  Dashed  lines  correspond  to  negative  values. 

Time  trace  along  centre-line  and  spanwlse  distribution 
at  intermediate  time  extracted  to  bottom  and  right  of 
plot  respectively. 

(a)  uncontrolled  reference  case 

(b)  c(xitrolled  case 
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Fig.  1 


EXTERNAL  CONTROL  -  THREE  FREQUENCIES 


r*.yi 


(b) 
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2  Centre-line  hot-wire  trace  at  target  location. 
Excitation  with  three  discrete  frequencies. 

(a)  uncontrolled  case 

(b)  "external"  control  on 
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Fig.  1 


SPANHISE  DISTRIBUTION 


*  Forcing  at  Re,.  ■  870,  control  at  Re^. 
target  zone  at  Re,,  «  1780 


1260, 


3  Integrated  hot-wire  ms  in  T-S  frequency  band  at  target 
location  versus  spanwise  probe  position. 

"External"  control  of  two  superposed  frequencies. 
Actuator  location  at  control  locaticxi  indicated  by 
full  triangles  on  abcissa. 


Fig.  14  Hot-wire  trace  at  detector  location  (centred  in 

atreamwise  position  between  double  microphone  but 
displaced  sideways  by  8mm  to  avoid  interference). 
Pulse  excitation  near  leading  edge. 

(a)  single  realization  hot-wire  record 

(b)  ensemble  average  of  64  realizations 
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Fig.  15  Microphone  trace  at  detector  locaticm  on  centre-line. 
Pulse  excitation  near  leading  edge.  Simultaneously 
recorded  with  hot-wire  trace  in  Fig.  14a. 

(a)  Single  realization  of  single  microphone 

(b)  Single  realization  of  difference  micrc^one 
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Fig.  16  Centre-line  hot-wire  trace  at  target  location. 
Single  mode  excitation  neau:  leading  edge. 

(a)  uncontrolled  signal 

(b)  controlled  signal  using  control  loop. 
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Fig.  17  Actuator /detector  cross-talk.  Detector  output  as 
a  function  of  externally  driven  actuator  forcing 
level . 

(a)  two  actuators  symmetrical  to  detector  microphones, 
spanwise  displaced  by  5mm. 

(b)  two  actuators  s^nmetrical  to  detector  microphones, 
spanwise  displaced  by  15mm. 
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EXPERIWENTAI.  SETUP 


*  One  and  twin  source  forcing. 

*  Measuring  station  at  1235. 

Cenfrot  tQr^<t  at  *  1780. 


Fig.  Cl  Aluminium  test  plate  for  active  control  experiments. 


EXPERIMENTAI.  FACILITY 


-  C.U.E.D.  low  turbulence  research  tunnel 
(former  NPL  tunnel). 

*  Exchangeable  working  sections  3ft  x  3ft  x  6ft. 

*  3-D  computer  controlled  traverse 

*  Micro-computer  based  d.a.  system 

*  Turbulence  Intensity  <  0.01% 

(20m/s,  4  Hz  to  2000  Hz). 


Fig.  C2  Cambridge  University  Low  Turbulence  Research  Tunnel. 
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Fig.  C3  Electronic  hardware  for  active  control  experiment. 


0 


(b) 


0.13  secs 


Fig.  D1  Boundary  layer  response  in  packet  centre-line  at  Re^  1380. 
Data  analogoue  filtered  between  20Hz  and  2000Hz. 

(a)  single  realization 

(b)  ensemble  average  of  189  realizations 
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Fig.  D2  Power  spectral  density  versus  frequency  plot  of  ensemble 
averaged  data  shown  in  Fig lb. 


1 


0 


25  msec 


Fig.  D4  Measured  u-fluctuaticx»s  in  the  boundary  layer  at  Re^.^  1235 
for  twin  pulse  interaction  experunent.  Two  pulses  of  same 
strength  buth  opposite  sign  excite  flow  at  Re^=  880. 
Ensemble  averaged  data  digitally  filtered  between  80Hz  and 
500Hz . 

(a)  perspective  view 

(b)  contour  plot,  negative  contours  dashed 


18.6  cm 


25  msec 


0  0 


18.6  cm 


25  msec 


Fig.  D5  Measured  u- fluctuations  in  the  boundary  layer  at  1235 
for  triple  pulse  interaction  experiment.  Three  pulses  of 
same  strength  (3cm  apart  spanwise)  excite  flow  at  Re^*=  880 
Middle  pulse  of  opposite  sign.  Ensemble  averaged  data 
digitally  filtered  between  80Hz  and  500Hz. 

(a)  perspective  view 

(b)  contour  plot,  negative  contours  dashed 
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Appendix  B:  Project  Plan 


Phase  Che  -  FYellniinaries 

-  Review  of  Literature. 

-  Familiarization  (R.H.)  with  wind  tunnel  setup,  computer 
controlled  data-acquisition  system. 

-  Test  plate  design. 

-  Three-dimensional  traverse. 


Phase  Two  -  Transfer  fincdons 

-  Selection  of  method  for  detection  and  control . 

-  Experimental  determination  of  detector  euid  controller 
tr2msfer  function. 

Phase  Three  -  Numerical  rrvddling 

-  Implement  2md  experiment  with  computer  model  of  various  control 
arrangements  using  Caster's  (1978)  rapid  series  eigenvalue 
expansion  method. 

-  Based  on  above  modelling  identify  the  most  promising  control 
arrangements . 

Phase  Four  -  Active  control  expaiment 

-  Experimental  study  of  the  most  promising  control  arrangements 
identified  in  phase  three  using  a  small  number  of  detectors 
and  controllers. 

-  Prepare  final  report. 
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Appendix  C:  Hardware  Development 

'nie  experiments  to  be  performed  require  a  flat  plate  which  is 
suited  for  the  testing  of  several  different  exciter,  detector  and 
controller  arrangements.  A  plate  (Fig.  C1 )  which  allows  for 
chemgeable  inserts  has  been  designed,  build  and  installed  in  the 
Cambridge  University  Engineering  Department  (C.U.E.D.)  low 
turbulence  wind  tunnel  (Fig.  C2).  The  plate  allows  disturbance 
excitation  at  five  spanwise  locati(xis  (3cm  apart)  via  embedded 
loudspeakers  20.5cm  downstream  of  the  leading  edge. 

An  importeint  feature  of  the  plate  is  the  provision  for  exch2uigeable 
inserts  which  allow  the  testing  of  various  different  control 
airrangements .  Inserts  can  be  exch^mged  while  the  plate  is 
positioned  in  the  tunnel.  Great  C2u:e  was  taken  that  the  edge  of  the 
insert  was  flush  with  the  working  side  of  the  plate  in  order  to 
avoid  the  creation  of  flow  disturbances.  The  gap  between  insert 
and  plate  was  filled  and  polished  in  order  to  achieve  the  required 
surface  continuity. 

Control  actuators  and  disturbance  detectors  had  to  be  selected 
according  to  frequency  response  characteristic,  ease  of 
installation,  physical  robustness  and  availability.  The  ’buried 
loudspeaker"  technique  had  been  successfully  used  in  previous  work 
(Caster  and  Grant,  1975)  for  disturbance  excitation  in  a  similar 
setting.  It  was  decided  to  use  the  saune  technique  for  control 
actuation.  In  order  to  allow  close  spanwise  spacing  of  the  control 
actuators  a  physically  small  speaker  was  required.  Various 
earphone  speakers  were  tested.  Ihe  mini-speakers  EH3030  of 
Knowles  Electronics  Co.  (73  Victoria  Road,  Burgess  Hill,  West 
Sussex,  RH15  9LP,  Englemd)  were  found  to  be  the  most  suitable  in 
view  of  their  frequency  response  ctuu'acteristics,  their  relative 
high  power  output  and  their  small  [^ysical  dimensicxis  (7mm  x  3mm  x 
3mm). 
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Disturbance  detection  was  originally  envisaged  to  be  done  with 
either  hot-fllm  sensors  flush  mounted  on  the  plate  surface  or 
embedded  microphones  a  la  Kendall  (1990).  Hie  hot-film  sensors, 
however,  were  eliminated  because  of  mecluuiical  difficulties.  It 
was  very  time  consuming  to  attach  flying  leads  to  the  sensor 
without  producing  protrusions  on  the  plate  working  side  which 
cause  unacceptable  flow  interferences.  In  addition,  the  electronic 
requirements,  i.e.  a  bridge  circuit  for  each  sensor  and  possibly 
drift  compensation,  is  rather  complex.  Hence,  embedded  microphone 
detectors  were  used.  The  micro^^xxie  detector  chosen  was  a  9.5mm 
diameter,  10mm  high  model  TANDY  No. 270-090.  The  detector 
communicated  with  the  plate  working  side  via  a  1mm  diameter  hole. 
Voltage  output  reading  was  readily  available.  The  output  signal 
level  was  only  about  0.2mV  amd  cau-eful  screening  was  required. 


Electronic  Hardware 

Figure  C3  shows  diagramatically  the  electronic  control  cirucuit 
employed  in  the  final  experiments.  Excluding  the  integrator  which 
was  only  used  during  part  of  the  experiment,  and  the  precision 
shielded  op-amps,  the  components  employed  were  standaurd  of  the 
shelf  compcxients.  The  time  constants  of  the  Bessel  Filters  were  to 
large  to  respond  to  a  passing  wave  packet  in  time  and  were 
partially  responsible  for  the  failure  to  achieve  wave  packet 
control.  The  Phase  Shifters  caused  narrow  bandwidth  filtering 
around  the  principle  frequency  leading  to  control  actuator  ringing 
after  the  detector  signal  had  ceased  to  exist. 
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Appendix  D:  Preliminary  Wave  Packet  Experiment 

Using  the  newly  designed  test  plate,  some  preliminary  wave  packet 
experiments  were  carried  out.  The  purpose  of  these  experiments  was 
threefold:  (a)  to  show  that  wave  packets  could  be  generated  with 
the  embedded  loudspeaker  technique,  (b)  to  demonstrate  that  the 
background  disturbance  level  in  the  C.U.E.D.  low  turbulence 
facility  was  indeed  low  enough  to  allow  easy  detection  of  a  single 
wave  packet  in  the  boundeury  layer  auid  (c)  to  investigate  if 
loudspeaker  excitation  was  "reversible,"  i.e.,  does  a  sign  chamge  in 
the  driving  voltage  pulse  create  a  wave  packet  of  same  amplitude 
but  180  out  of  phase?  Point  (b)  is  especially  important  for  the 
implemer  ^ation  of  the  control  system  since  some  filtering  will 
have  to  be  employed  to  isolate  the  relevauit  wave  packet  features 
from  "noise"  present  in  the  boundary  layer.  In  addition  to  these 
three  main  objectives,  the  experiments  were  required  to  test  the 
new  plate  setup  for  (near)  zero  streamwise  pressure  gradient  and  to 
check  the  three-dimensional  traverse  gear  and  data  acquisition 
software. 

Two  different  sets  of  experiments  were  performed.  The  tunnel  speed 
for  all  experiments  was  20m/sec  amd  measurements  were  taken  using 
a  constant  temperature  boundary  layer  hot-wire.  A  "near  zero" 
streamwise  pressure  gradient  was  achieved  by  adjusting  a  flap 
extending  downstream  from  the  main  plate.  In  order  to  avoid  flow 
separation  at  the  leading  edge  a  second  small  trailing  edge  flap 
was  adjusted  to  move  the  stagnation  point  slightly  to  the  working 


-  67  - 


side  of  the  plate.  Following  work  by  Caster  and  Grant  (1975),  the 
hot-wire  was  positioned  just  outside  the  boundary  layer  at  1 .2 
boundary  layer  thicknesses  from  the  plate  surface,  where  the 
Tollmien-Schlichting  eigenfunctions  show  an  outer  (flat)  meucimum. 

In  the  first  set  of  experiments,  the  downstreaun  and  spanwise 
development  of  a  single  wave  packet  was  investigated.  The  exciter 
location  at  20. Sera  from  the  leading  edge  corresponds  to  Re  =  880 
(Reynolds  number  based  on  displacement  thic)cness)  and  the  flow 
responses  were  measured  at  R^=  1235  (40.5cm)  and  1380  (50.5cm), 
respectively.  Typically  45  spanwise  "data-slices"  roughly  six 
displacement  thicknesses  apart  were  recorded.  The  exciting  pulse 
was  computer  generated  and  could  be  varied  in  cunplitude  and  width. 
For  the  experiments  described  in  this  report  a  constant  pulse  width 
of  0.5msec  was  used.  To  allow  ensemble  averaging  of  the  data,  the 
pulse  generation  was  synchronised  with  the  data  acejuisition  system. 
Itie  sampling  rate  was  4kHz  and  the  signal  was  analogue  filtered 
between  20Hz  2md  2kHz.  Each  record  contains  512  data  points  and 
collection  time  per  record  was  thus  0.13sec.  An  ensemble  average 
was  formed  using  189  individual  records. 

Figure  Dia  shows  a  typical  single  realization  of  the  boundary  layer 
response  at  Re^fs  1380.  Figure  Dib  shows  the  same  data  after 
ensemble  averaging.  A  power  spectral  density  versus  frequency  plot 
of  the  ensemble  averaged  data  is  shown  in  Fig.  D2.  The  quite 
satisfactory  signal -to- noise  ratio  is  apparent.  Ihe  Tollmien- 
Schlichting  frequency  band  shows  power  spectral  densities  of  about 
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40dB  above  the  background  noise  level.  A  perspective  view  and  a 
contour  plot  of  the  excited  wave  packet  at  the  two  different 
Reynolds  numbers  are  shown  in  Figs.  4  and  D3.  In  order  to  emphasize 
the  Tollmien-Schlichting  wave  response,  the  ensemble  averaged 
records  presented  in  the  last  two  figures  were  banned  and  digitally 
filtered  between  80Hz  2md  500Hz.  At  the  higher  Reynolds  number, 
the  stronger  three-dimensional  modulation  of  the  signal  aund  the 
increased  packet  "length"  are  apparent. 

The  first  set  of  experiments  established  that  (a)  the  embedded 
loudspeaker  technique  is  suited  to  the  generaticm  of  wave  packets 
in  laminar  boundary  layers  and  (b)  the  C.U.E.D.  low  turbulence  wind 
tunnel  provides  the  desired  flow  quality  for  this  type  of 
experiments. 

"rtie  second  series  of  experiments  involved  the  interaction  of  wave 
packets  generated  simultaneously  at  different  spanwise  locations. 
Only  the  signal  at  Re^.*  1235  was  recorded.  First,  two  pulses  3cm 
apart  of  same  amplitude  but  opposite  sign  were  introduced  at  the 
exciter  location  to  disturb  the  flow,  llie  experiment  was  then 
repeated  using  three  spanwise  drivers  3cm  apart.  Again  all  pulses 
generated  had  the  same  amplitude  but  the  middle  pulse  was  of 
opposite  sign.  The  data  acquisition  and  data  reduction  performed 
was  as  described  above.  Here,  however,  55  and  63  spanwise  "slices" 
were  collected  in  order  to  map  out  the  interaction  region  for  the 
"twin"  pulse  and  "triple"  pulse  experiments  respectively. 
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Figures  D4  and  D5  show  perspective  views  and  contour  plots  of  the 
flow  response  to  twin  pulse  and  triple  pulse  excitation.  Itie 
quality  of  the  data  obtained  is  not  as  high  as  in  the  first  set  of 
experiments  reported  above.  This  was  especially  so  in  the  triple 
pulse  experiment  which  showed  2ui  underlying  modulation  of  about 
500Hz.  This  was  traced  to  vibrations  of  the  hot-wire  when  placed  in 
a  shearflow,  a  phenomenon  traced  to  lack  of  wire  tension.  Ibe 
contour  plots  of  both  experiments  show  the  symmetry  of  the  flow 
response  for  "positive"  and  "negative"  pulse  excitation.  It  should 
be  pointed  out  that  the  lines  represent  the  seune  absolute  values 
for  positive  zmd  negative  contours  drawn.  The  contour  plots  for 
both  experiments  show  nairrow  regions  of  wave  cancellation  between 
adjacent  wave  packets  of  opposite  sign. 

The  second  series  of  experiments  achieved  objective  (c)  formulated 
above,  i.e.,  sign  switching  of  the  driving  pulse  leads  to  sign 
switching  of  the  resulting  wave  packet. 
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Appendix  E:  Programme  Listings 
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Apf^cudtK  £•  UtiiM: 

Prograun  maincon  *4 

************************************************************ 

Control  of  instability  waves 

march  31,  91  modified  to  allow  phase  shifting  in 

order  to  asses  sensitivity  to  off-timing 

changed  to  read  in  unformatted  data  files  to  speed  up 

run  (march  29,  15:00) 

july  16,  91  maincontr2.f  with  double  lo<^  (110,105) 

to  map  out  ampgain,  phaseshift  space  (Hch) . 

jxaly  18,  91  maincontr3.f  implements  one  detector  driving 
several  (three)  actuators  spemwise  where  the 
two  sideactuators  are  driven  with  anplitude 
weight  of  0.5  of  the  middle  actuator  amplitude  (Hch) 

july  25,  91  maincontr4.f  ia^lements  random  forcing  in 

space  and  time.  Note  that  FUNCTAION  DRAN2  had 
to  be  corrected. (Hch) 

july  26,  91  maincontrS.f  speeds  v?)  convolution  with 
known  system  tramsfer  functions  by  using 
option  of  reading  unformatted  Fourier 
coefficients  of  these  from  data  directory 
DAT/Fsy....  Hence  convol2d  is  modified 
in  order  to  take  advantage  of  this.  (Hch) 

Furthermore  pogram  was  modulized  into 
more  subroutines  and  common  block 
is  used  to  reduced  compiled  file  storage 
requirements. 

sept  5,  91  modified  just  for  target  field  superposition 
of  vhole  fields  (in  scontr.f)  (maincontr6.f) 

Subroutines  required: 

sin  -  prepare  ^ysical  space  data  field 

at  location  of  microphone  detectors 
written  in  array  Hei^tl 

convol2d  -  2D  feist  Fourier  transform  of  512x128 

double  precision  fields.  If  both  fields 
aire  in  physical  space  (ini  =  in2  =  0),  Hei^tl 
then  physical  space  result  is  in  Heightl. 

If  arrays  are  in  proper  fourier  form  (ini  and/ 
in2  -  1)  then  array  tempi  and/or  teiip2  are  use 
and  convolution  is  speeded  up  significantly. 
Convoluted  result  is  always  returned  in  Height 
(physical  space) . 

FOURN  -  ID  fast  Fourier  transform 

scontr  -  control  of  distuzt>anoe  field  and  continuaticxi 
of  controlled  field  to  teuiget  location. 
Uncontrolled  physical  space  input  field  in 
Heightl,  controlled  ouput  field  also  in 
Heightl  (Height2  is  dimay  field  -destroyed) . 

sanaly  -  Analyze  control  performance  at  target  location 
in  certain  unspoiled  part  (wraparound)  of 
field  Height2. 

sout  -  output  of  resulting  field  to  haxddisc  director 
. /DAT/<f ilename> 

Data  files  from  directory  ./DAT: 

Fsystrl  -  file  with  Fourier  space  representation 

of  system  tremsfer  function  I  (unformatted 
double  precision) . 
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Fsystrll  -  file  with  Fourier  space  representation 

of  system  transfer  function  II  (double  precis! 
unformatted) 


Roland  Heinrich 


Ceunbridge  , March  27,  1991 


c  ************************************************************* 

c 

implicit  real{a-b,d-h,o-y)  ,conplex(c,z)  ,integer(i-n) 

Character*  16  fnaune 
Character*!  A 

c  parameter  (inmx=10,iampmx=2  2) 

parameter  ( inmx=l ,  iampmx=2 ) 

parameter  (ndim  =  2,itmzoc=512,k2roax=128,ivec=itmax*kzmax*2) 
double  precision  psen,ai!pgain,pi,  factl,  fact2,  fact3,  fscale,help, 
+  ampdel 

dimension  psvec(l:  ieunpmx, -inmx:  inmx) 
c 

double  precision  tempi, tenp2,Heightl,Height2 
common/fld/teiipl  ( ivec) ,  temp2  (ivec) , 

+  Heightl{itmax,kzm^oc)  ,Height2  (itmax,kzmax) 

c 

pi  =  4.dO*datan(l.dO) 

print*, 'pi  =  ',pi 

idum  =  -1 

fname(l;4)  ='QAT/' 
c 

c  *****  set  phys.  space  disturbcince  field  at  Detector  location 
c 

c  1  call  sinf 

c 

c  *****  control  by  setting  selective  spanwise  values  to  zero 
c  centerline  is  at  l,use  1  <  ncontrol  <  64 

c 

1  print*,'  What  control  amplifier  gain  ?  ' 

read(*,*)air.pgain 

print*,'  Timeshift  nshift*.25msecs,  nshift  =  ?' 
read(*,  *)n£>hift 


print*,'- - ' 

print* ,  'au  pgain  ' ,  eunpgain 
nshift  =  i  s 

print*, 'ns ni ft  ', nshift 
call  sconti  (ampgain, nshift) 
goto  555 
c 

c  *****  Pull  whole  field  forward  by  npull  such  that  noncorrupted  field 
c  (by  wraparound)  starts  at  ntime  =1  to  ntime  =  420. 

c  Note  also  that  in  spanwise  direction  noncorrupted 

c  field  lies  between  iz=35  2Uk1  iz=93.  Importzmt  for  pcwer 

c  calculation. 


npul  =  360 

do  ik  =  l,kzmax 

ii  =0 

do  it  =  npul,itmzo( 
ii  =  ii  +  1 

Heic^t2(ii,ik)  =  Heightl(it,ik) 
end  do 
end  do 

do  ik  =  l,kzroax 

ii  =  itmax-rpul+1 

do  it  =  l,rpul-l 

ii  =  ii+1 

Height2(ii,ik)  =  Heightl(it,  ik) 
end  do 
end  do 


c  *****  tinnlysr>  rnntro]  porformancp  it  tar-ir-t 


in  un-spoiled  area  of  field  Height2 


c 
c 

itl  =  1 
ith  =  420 
izl  =  35 
izh  =  93 
555  itl  «  1 

ith  «  512 
izl  =  1 
izh  *=  128 
do  ik  =  Ijkzmax 
do  it  =  Ijitmax 

Height2  (it,  ik)  =  Heightl(it,ik) 
end  do 
end  do 
c 

call  sanaly(itl,ith,izl,izh) 
c 

c  *****  write  out  result 
c 

call  sout (Height2) 
c 

pause  'new  control  configxiration?' 
goto  1 

900  print*, 'Error  in  reading  ',fncune,'  start  over  ?' 
pause 
goto  1 

910  print*, 'End  reached  in  reading  ',fnaine,'  start  over?' 

pause 
gotol 

999  end 

Subroutine  sinf 

c  ************************************************************** 

c 

c  prepares  Heightl(itmax,kzniax)  to  correspond  to  the  uncontrolled 
c  disturtjance  field  at  the  location  of  the  detector  (and 
c  perharps  controller  array) . 

c  Outpjut  :  Heightl{itBax,kzinax) 

c 

c  Subroutines  required;  convol2d 

c  Foum 

c 

c  Roleuid  Heinrich  Ceunbridge,  July  31,  91 

c  ************************************************************** 

Character*16  fneune 
Character*!  A 

parameter  (itineoc=512  ,kzincix=128 ,  ivec=itm20c*kzmax*2) 
double  precision  tempi, temp2,Heightl,Hei^t2 
ccraiDon/f Id/ tempi  ( ivec ) ,  tenp2  ( ivec) , 

+  Heightl(itmax,kzmax)  ,Height2(itinax,kzmax) 

c 

fname(l:4)  =  'DAT/' 

1  print*,'  Do  you  have  double  precision  phys.  space  ' 

print*,'  input  file  corresponding  to  disturbcuxre  ' 
print*,'  flow  field  at  control  location  (1)  or  do  ' 
print*,'  you  have  Fourier  space  input  data  (2)  at  ' 
print*,'  exciter  location  (all  product  of  ranfield' 
print*,'  input  1  or  2  :  ?' 
rcad(*,*)  infld 
if (infld.eq.l)  then 

print*,'  Fhys.  dist.  field  at  control  location  (doubl.prec.)? 
rcad(*,*)  fna]ne(5:16) 
print*, f name 

open  ( 1 ,  f  ile=fname,  forro="unforraatted" ) 
read(l,ERR=900,ENI>=910)  Heightl 
close (1) 
goto  10 

print*, 'Error  in  reading  ',fnaine,'  try  again  ' 

pause 


900 
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goto  1 

910  print*, 'End  of  file  encountered  in  ',fneune,'  try  again' 
pause 
goto  1 

else  if  (infld.eq.2)  then 

open(l,  file="DAT/FsystrIun*' ,  forro="unformatted") 
read(l,ERR=900,E3»I>=910)  tenp2 
close (1) 
c 

c  ****  Read  in  fourier  trzmsform  of  input  field  (512x128) 
c  which  was  generated  with  rauifield.f 

c 

print*, '  Fourier  space  field  of  disturbemce  input  ?  ' 
read(*,*)  fneuDe(5:16) 
open(l,  file=fna]iie,  forn»=“unfonnatted'') 
read(l)  tenpl 
close (1) 
c 

c  *****  convolute  zcmdom  field  with  systrl 
c 

inl=l 

in2=l 

call  convol2d(inl,in2) 
else 

pause  'Wrong  infield  integer,  try  again' 

goto  1 

endif 

10  print*, 'Phys.  dist.  field  at  detector  location' 
print*, ' (Heightl)  set. ' 
c 

return 

end 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


subroutine  convol2d  ( ini ,  in2 ) 

************************************************************* 

Use  2-d  Fast  Fourier  transform  (PCXJRN)  to  convolute  arrays 
of  itni2oc  t-points  *  kzin20c  z-point  of  data. 

July  26,  1991  modified  to  allow  convolution  and  backtransform 
if  fields  are  passed  dcwn  as  appropriate  fourier 
coefficients  in  vector  tenpl  and/or  teiip2  eis 
indicated  by  control  Iparameter  ini  and  in2. 
ini  =  0  :  data  in  datinl  need  transform 
ini  =  1  :  fourier  tranformed  data  in  tenpl 
are  ready  for  multiplication. 
in2  =  0  euid  in2  «  1  as  above  for  ini.  (Hch) 

tenpl, datinl  :  holding  convoluted  result 
upon  return 

temp2,datin2  :  unaltered  by  subroutine 

Rolemd  Heinrich  Cambridge,  march  26,91 

**************************************************** 

implicit  real  (a-h,o-y) 
double  precision  pi,helpr,helpi 
character*12  fname 

parameter (ndiffl  =  2,itmax  -512,kzm20c=128,ivec=itmax*)czmax*2} 
dimension  nvec(ndim) 

double  precision  tenpl, tenp2, datinl, datin2 
common/ fid/ tempi  ( ivec) ,  tenp2  ( ivec) , 

-t-  datinl  ( itmax,  kzmax)  ,datin2(itinax,)czmax) 

nvec(l)  =  itmax 
nvec(2)  =  kznmuc 


c 


pi4  =  datan(-l.dO) 


o  o 


pi  =  4. do*  pi4 
c 

c  ***  prepare  complex  data  vectors  for  subroutine  POURN 
c  if  necessary  as  indicated  by  dummy  argument  inl,in2 
c 

If(inl.eq.O)  then 
jj  =  1 

do  10  ik  =  l,k2oax 
do  15  it  =  ijitmax 

tenpl(jj)  =  datinl(it,ik) 
templ(jj+l)=  O.dO 

jj  =  jj  +  2 

15  continue 

10  continue 
c 

c*  check  nax  dimension 

if  ( (jj-1)  .ne.ivec)  stop  'convol2d  jj  n.e.  ivec  ' 
c 

c  ***  transform  into  fourier 

call  FOURN(teBpl,nvec,ndim,+l) 
end  if 
c 

If(in2.eq.O)  then 
jj  =  1 

do  11  ik  =  l,kzinax 
do  17  it  =  l,itmax 

teBp2(jj)  =  datin2  (it,ik) 
teiiip2(jj+l)=  O.dO 

jj  =  jj  +  2 

17  continue 

11  continue 
c 

c*  check  max  dimension 

if  ( (jj-1)  .ne.ivec)  stop  'convol2d  jj  n.e.  ivec  ' 

***  transform  into  fourier 

call  P0URN(tenp2,nvec,ndim,+l) 
end  if 
c 

c  ***  convolute  data  in  fourier  space,  teii5)l(jj)  holds  convolution 
c  this  is  start  point  if  ini  and  in2  .ne.  0  indicating  that 
c  tempi, temp2  hold  allready  fourier  coefficients  in  proper  form 
c 

do  16  jj  =  1, ivec, 2 

helpr  =  tempi (jj)*temp2(jj)  -  tempi (jj+1) *temp2(jj+l) 
helpi  =  teii?>l(jj)*ten?}2(jj+l)  +  ten?)2(jj)*teii5>l(jj+l) 
templ(jj)  =  helpr 
tempi  (jj+1)  =  helpi 

c  print*, 'convoluted,  j,j+l,',  jj,templ(jj)  ,ten?)l(jj+l) 

16  continue 
c 

c  ***  backtransform  into  physical  spac^ 
c 

call  FOURN (tempi, nvec,ndim,-l) 
fact  =  l.d0/dfloat(nvec(l)*nvec(2) ) 
c 

c  *****  rearrange  data  vector  in  (time, space)  array 
c  (see  arreingement  convention  'Num.  Recipes'  pg  451) 

c  Note  that  the  imaginary  part  of  convoluted  sum  in 

c  physical  space  should  be  identical  zero.  Check! 

c 

jj  =  1 

dsumim  =  O.dO 
do  20  ik  =  l,kzmax 
do  25  it  =  l,itroax 
datinl(it,ik)  =  fact*templ(  jj) 
dsumim  =  dsumim  +  dad)s(templ(j  j+1) ) 

jj  =  jj  +  2 

25  continue 
20  continue 
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print*,'  SUB  of  imag.  part  should  be  zero  ',dsuttiB 
return 
end 


Subroutine  FXXJRN(data,nn,ndiin,isign) 

c**  ***************************************************** 

c 

c  Replaces  DATA  by  its  NDIM-dimensional  discrete  Fourier 
c  tremsform  if  ISIGN  is  input  as  1.  NN  is  an  integer 
c  array  of  length  NDIM,  containing  the  lengths  of  each  dimension 
c  (number  of  conplex  values) ,  which  MUST  all  be  powers  of  2. 
c  DATA  is  a  real  array  of  length  twice  the  product  of  these 

c  lengths,  in  which  the  data  are  stored  as  in  a  multidimen- 

c  sional  conplex  FORTRAN  array.  I  ISIGN  is  input  as  -1,  DATA 

c  is  replaced  by  its  inverse  transform  times  the  product  of 

c  the  lengths  of  all  dimension, 
c 

c  from  Press,  Flannery,  Teu)colsky,  amd  Vetterling, 
c  "Numerical  Recipes",  Ccunbridge  University  Press,  !986, 

c  pg.  451. 

c 

c  Converted  to  overall  double  precision  (R. Heinrich) 
c 

c  Roleuid  Heinrich  Caunbridge,  Febr.  8,  1991 

c******************************************************* 


c 

Implicit  real  (a,b,d-h,o-y) ,  complex(c,z) , 

+  integer (i-n) 

Double  Precision  wr,wi,wpr,wpi,wtenp, the ta,tempr, tempi, 
+  data 

Dimension  nn(ndim)  ,data(*) 
c 

c  ***  Compute  total  number  of  complex  values, 
c 

ntot  =  1 

do  11  Idim  =  l,ndim 
ntot  =  ntot*nn(idim) 

11  continue 

print*, 'foum  ntot,ndim  ',ntot,ndim 
c 

c  do  500  i  =  l,ntot*2 

c  print*, 'foum-i, data  ',i,data(i) 

c  500  continue 
c 

c  ***  Main  loop  over  the  dimensions, 
c 


c 

c 

c 


nprev  =  1 

do  18  idim  =  l,ndim 
n  =  nn(idim) 
nrem=  ntot/ (n*nprev) 
ipl  =  2*nprev 
ip2  =  ipl*n 
ip3  =  ip2*nrem 

***  This  is  the  bit  reversal  section  of  the  routine 


12 


i2rev  *=  1 

do  14  i2»l,ip2,ipl 
if (i2.lt. i2rev)  then 
do  13  il=i2,i2+ipl 
do  12  i3=il,ip3, 
i3rev 
tenqpr 
temoi 
data(i3) 
data(i3-»-l) 
data(i3rev) 
data(i3rev+l) 
continue 


-2,2 

ip2 

*  i2rev+i3-i2 
=  data(i3) 

=  data(i3+l) 

=  data(i3rev) 

=  dataU3rev+l) 
=  tempr 
=  tempi 
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l 


13  cxjntinue 
endif 

c 

ibit  =  ip2/2 

1  if  ((ibit.ge.ipl)  .and.  (i2rev.gt.ibit))  then 

i2rev  =  i2rev-ibit 
ibit  =  ibit/2 
goto  1 
endif 

i2rev  =  i2re!v  +  ibit 

14  continue 
c 

c  ***  Here  begins  the  Danielson-Lanczos  section  of 
c  the  routine . 
c 

ifpl  =  ipl 

2  if (ifpl.lt.ip2)  then 

ifp2  =  2  *  ifpl 
c 

c  ***  Initialize  for  trig  recurrence 
c 

theta  =  isign  *  6.28318530717959d0/ (ifp2/ipl) 
wpr  =  -2.d0*dsin(0.5cj0*theta)  **2 

wpi  =  dsin (theta) 

wr  =  l.do 

wi  =  O.dO 

c 

do  17  i3  =  1, ifpl, ipl 
do  16  il  =  i3,i3+ipl“2,2 
do  15  i2  =  il,ip3,ifp2 
kl 
k2 

tenpr  « 

tempi  ® 

data(k2)  = 
data(k2+l)  = 
data(kl)  == 
data(kl+l)  * 

15  continue 

16  continue 
c 

c  ***  trigonometric  recurrence, 
c 

wtenp  =  wr 
wr  =  wr*wpr  - 

wi  =  wi*wpr  + 

17  continue 
c 

ifpl  =  ifp2 
goto  2 
endif 

r^rev  =  n*nprev 

18  continue 
return 

end 

Subroutine  scontr(ampgain,ni) 

c  **************************************************************** 

c  Controls  field  according  to  user  interaction  with  regard  to 

c  control  inputs  and  control  anqplitude  and  possible  phase  shift 

c  fieldl  =  uncontrolled  input  field 

c  holds  controlled  field  iqxm  return 

c  field2  -  dummy  field  which  is  destroyed  in  routine 

c 

c  Roland  Heinrich  Cambridge,  july  29,  1991 

c  **************************************************************** 

c 

character*16  fname 

double  precision  ampgain,help 

parameterfndim  =  7,  itmax=512  ,  kzmay-  1  28  ,  ivr>c 


i2 

kl  +  ifpl 

wr*data  (k2 )  -wi*data  {k2+l) 
wr*d4ta (k2+l ) +wi*data (k2 ) 
data(kl)~tempr 
data(kl+l)  -  tempi 
data(kl)  +  tempr 
data(kl+l)  +  tempi 


wi*wpi  +  wr 
wteiqp*wpi  +  wi 
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c 

double  precision  ten?)!, tenp2,Fieldl,Field2 
common/ fid/ teapl  ( ivec) ,  teinp2  ( ivec) , 

+  Fieldl  ( itmax ,  kzmax) ,  Field2  ( itmax ,  kzmax) 

c 

c  *****  preset  Field2  with  O.dO 
c 

do  ik  =  1, kzmax 
do  it  ■=  l,itmeoc 
Field2(it,ik)  =  O.dO 
end  do 
end  do 
c 

c  *****  control  by  setting  selective  spcinwise  values  to  zero 
c  centerline  is  at  64, use  1  <  ncontrol  <  128 

c  Control  actuator  input  is  2*3inm  to  either  side 

c  of  detector  line  with  wei^t  0.5 

c 

print*,  ^Control  phase  -  :  Set  selected  microphone  position' 
print*,'  Actuate  at  position  +/-  6nim.' 

90  print*, 'Enter  micontrol  >  1000  if  you' 

print*, 'have  set  all  control  points.  3  <=  ncontrol  <=126. ' 
print*,  'Enter  mic-control  =?  ' 
read* , ncontrol 
if  (ncontrol.eq.o)  then 
c  *****  no  control 

do  ik=  l,kzmcix 
do  it  =  1, itmax 
Field2(it,ik)  =  Fieldl ( it, ik) 
end  do 
end  do 
end  if 
c 

if  (ncontrol. GT.  1000)  goto  95 
c 

do  it  =  1, itmax 

help  =  Fieldl (it, ncontrol) 

Field2(it,ncontrol-2)  =  Field2( it, ncontrol -2)  -  0.5d0*help 

Field2(it,ncontrol+2)  =  Field2  (it,ncontrol+2)  -  0.5d0*help 

end  do 
c 

goto  90 
c 

c  convolute  controlled  Field2  with  Systrll  (Fieldl) 

c 

95  open(3,  file="DAT/FsystrIIun",form=”unformatted") 
read (3)  tempi 
close (3) 
c 

ini  =  1 
in2  =  0 

call  convol2d(inl, in2) 
c 

c  *****  amplifier  gain  and  off-timing  of  control  signal 
c  in  ni*0.25  msec 

c 

if(ni.gt.O)  then 
do  it  ■  l,itmax-ni 
do  ik  =  l.kzmauc 

Field2(it,ik)  =  Fieldl  (it+ni,ik) 
end  do 
end  do 
do  it  =  l,ni 

do  ik  =  l,kzm20c 

Field2  ( itroax-ni+it,  ik)  =  Fieldl  (it,  ik) 
end  do 
end  do 
c 

elseif  (ni.lt.O)  then 

c 
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do  it  “  1,-ni 

do  iX  «  l,kzln^lx 

Field2  (it,  ik)  =  Fieldl(itmao{+ni+it,ik) 
end  do 
end  do 
c 

do  it  =-ni+l, itmaoc 
do  ik  =  l,k2inax 

Field2(it,ik)  =  Fieldl(ni+it,  ik) 
end  do 
end  do 
c 

else 

c 

do  it  =  Ijitmax 
do  ik  =  l,kzinax 

Field2  (it,  ik)  =  Fieldl(it,ik) 
end  do 
end  do 
c 

end  if 
c 

c  *****  now  read  in  uncontrolled  field  \diich  is  to  be 
c  superposed  with  control  1  output 

c 

200  fnetine='DAT/ptargrcin' 

open  ( 3 ,  f  ile=f  name ,  f  onn=”unf  ormatted" ) 
read(3,ERR=900,END=910)  Fieldl 
close (3) 
c 

196  do  ik  =  l,kzmax 

do  it  =  l,itmax 

Fieldl  (it,  ik)  =  ampgain*Field2(it,ik) 

1  +  Fieldl (it, ik) 

end  do 
end  do 
return 

900  print*, 'ERRor  reading' ,fname, '  try  again' 
pause 
goto  200 

910  print*, 'END  of  file  ',fncune,'  try  again  ?' 
pause 
goto  200 
end 


Subroutine  s€maly(itl,ith,izl,izh) 

c  **************************************************************** 

c  Analyze  data  in  in  subset  itl-ith,izl-izh  of  field2  itmax  x  kzin20c 
c  (  pseudopower  per  point  ) 

c  field2  =  controlled  input  field 

c 

c  Rolemd  Heinrich  Cambridge,  july  30,  1991 

c  **************************************************************** 

c 

double  precision  psen,rpsen 

parameter  (ndim  =  2,itmax=512,kzmax=128,  ivec=iti»ax*kzmax*2) 
c 

double  precision  tempi, tenp2, Fieldl, Field2 
common/ fld/templ(ivec)  ,tenp2(ivec) , 

+  Fieldl ( itmax, kzmax)  ,Field2( itmax, kzmax) 

c 


c 

c 

c 


***  find  pseudo  disturbance  energy  in  field 


psen  =  O.dO 
icount  =  0 
do  ik  =  izl,izh 
do  it  =  itl,ith 
icount  =  icount+1 


.dc*; icidi 


end  do 

psen  =  dsqrt(psen) 
c 

c  *****  calculate  residual  power  per  point 

V. 

rpsen  =  psen/dfloat (icount) 

Print*,'  Power  in  controlled  subfield  ',psen 
Print*,'  Power  per  pcint  in  contr. field' ,  rpsen 
c 

cc  print*,'  psvec(ieunp,ni)  '  ,psvec(i2uip,ni) 

c 

ctest  go  to  101 

cl  110  continue 

cl  105  continue 

cll  print*,  '  Enter  filename  for  Power  vector  psvec (icunp,  in) 

c 

cl  fnaine(l:4)  =  'DAT/ ' 

cl  read(*,107)  fnaine(5:16) 

107  format (A12) 

cl  open  (2,file=fncUQe,fonD='formatted') 

cl  write(2,*)  ( (lamp, in, psvec (lamp, in) , 

cl  1  in=-innix ,  inmx) ,  iamp=l ,  ianpmx) 

102  close(2) 
c 

return 

end 

Subroutine  sout (Height  1) 

c  ************************************************************ 

c  Write  out  data  fields  for  possible  post-processing 

c  via  graphics  routines.  Hence  charnge  field  data 

c  to  single  precision  for  coun tourplotting  contsim*.f 

c 

c  Roland  Heinrich  Cambridge  , March  31  1991 

c 

C 

inplicit  real  (a-b,d-h,o-y)  ,coiiplex(c,z) ,  integer (i-n) 
Character*16  fname 
Character*!  A 

parameter (ndim  =  2,  itiiiax=512,kzmax=128,  ivec=itjnax*kzmax*2) 
double  precision  psen, psvec, Heightl 
dimension  Heightl (itmax,kznauc) 
dimension  sHeight (itmax,kzmax) 
c 

c  *****  convert  Height2  to  single  precision  for  later  use  in  contour 
c  plot  routine 

c 

do  ik  =  l,kzmeix 
do  it  =  l,itmeix 

sHeight  (it,  ik)  =  sngl  (Heightl  (it,  ik) ) 
end  do 
end  do 
c 

100  print*,'  Do  you  want  to  store  field  single  precision  (1),  ' 
print*,'  double  precision  (2)  or  both  (3)?  Enter  1,2  or  3; 
read(*,*)  istore 
if  ( istore. eq.l)  then 

print*,'  filename  to  write  controlled  result  (sngl.prec. ) ' 
print*,'  at  target  location  =  (in  subdir  DAT)' 
print*,'  (unformatted  store)  ' 
fname(l;4)  =  'DAT/' 
read* , fname (5:16) 

open (4 ,  f  ile=fneune,  ERR=900,  form="un formatted") 
write(4,ERR=900)  sHeight 
write(4,ERR=900)  psen 
close (4) 

else  if ( istore. eq. 2)  then 

print*,'  filename  to  write  controlled  result  (double  proc.) 
trint*,'  ;'t  t'nr'U''’''  ’rK',Tti'^ 
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print*,'  (unformatted  store)  ' 
fna]ne(l:4)  =  'DAT/' 
read* ,  fname (5:16) 

cpen  ( 4 ,  f  ile=f  neune ,  ERR=900 ,  f  orm="unf  ormatted" ) 
write(4,ERR=900)  Heightl 
write(4,ERR=900)  psen 
close (4) 

else  if  (istore.eq. 3)  then 

print*/  filencune  to  write  controlled  result  (sngl  .prec. ) ' 
print*/  at  target  location  =  (in  subdir  DAT)' 
print*/  (unformatted  store)  ' 
fneune(l:4)  =  'DAT/' 
read* , fname (5:16) 

open  ( 4 ,  f  ile=f  name,  ERR=900 ,  f  orm=”unformatted" ) 
write(4,ERR=900)  sHeight 
write(4,ERR=900)  psen 
close (4) 

print*/  filename  to  write  controlled  result  (double  prec.)' 

print* , '  at  target  location  =  (in  subdir  DAT) ' 

print*,'  (unformatted  store)  ' 

fname(l:4)  =  'DAT/' 

read* , fname (5:16) 

open  ( 4 ,  f  ile=f  neime ,  ERR=900 ,  form=”unf  ormatted" ) 

write(4,ERR=900)  Heightl 

write(4,ERR=900)  psen 

close (4) 

else 

pause  '  error  in  storage  integer  istore  ' 
goto  100 
end  if 

return 

print* /Error  in  write  operation  file  ', fname 

print*/  Try  again  ?  ' 

pause 

goto  100 

end 


900 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

1 


c 

c 


c 

c 


c 


program  remfld 

************************************************************ 


Prograimn  to  set  the  ii^t  data  field  (it]nax*kzmax) 
used  as  input  for  control  modelling  (maincontrS.f) . 

CXitput  is  written  out  either  as  physical  space  field 
at  detector  location  (this  requires  convolution)  with 
systrl  or  as  field  Fourier  component  at  exciter 
location.  In  either  case  is  the  resulting  file 
written  as  unformatted  double  precision. 

August  2,  1991  random  field  is  set  in  Fourier  space 
with  reindom  phase  2md  imit  amplitxide  of 
fourier  coefficients.  (Hch) 

ireuid  =  starting  integer  for  ramdcan  ntmaber  gen. 

Subroutines  required  :  FOURN 

sremfld 

dr2in2 

convol2d 

Roland  Heinrich  Cambridge  ,July  27  1991 

'klfk-k'kitifk-k'kifkitiK-k'kit'kifkitit'k-kic-k-k'k'kifklfk'k'kitifk'kifkiticit'kitifkltifkltitititititicieifie 

implicit  real  (a-b,d-h,o-y)  ,coi!55lex(c, z) , integer (i-n) 
Character*16  fnaime 
Character*!  A 

fiarameter  (ndinF=2 ,  itinax=512 ,  k2max=128 ,  ivec=itmax*kzmax*2 ) 
double  precision  tempi, temp2,datinl,datin2 
common/  fid/ tempi  ( ivec) ,  temp2  ( ivec) , 

+  datinl(itmax,>czmax)  ,datin2(itm2ix,)czmax) 

dimension  nvec(ndim) 

nvec(l)  =  itmax 

nvec(2)  =  kzmax 

irand  =  -1 

fname(l:4)='DAT/' 

C2dl  sreuifld( irand, tempi) 

print*,'  Enter  1  to  store  Fourier  coefficients  of  field' 
print*,'  at  exciter  location  or  2  to  calculate  phys.  field' 
print*,'  at  detoctor  location  downstream  ' 
read(*,*)  ifield 

if  (ifield. eq.l)  then 

print*, '*UN* formatted  file  for  Fourier  coef . ' 
print*,'  of  random  field  (double  prec.  )in  DAT  ?' 
read*,  fname(5:16) 
print* ,  f  neune 

open(l,file=fncune,  form="unfoirmatted") 
write (1)  tempi 
close (1) 

else  if  (ifield. eq. 2)  then 

open  (1 ,  f  ile="DAT/FsystrIun" ,  form=’'unformatted'') 
read(l)  temp2 
close (1) 

ini  =  1 
in2  =  1 

call  convol2d(inl, in2) 

print*, '*UN*forTOatted  file  for  phys.  field  at  ' 

print*,'  detector  location  downstream  (doubl.  prec.)  in  DAT?' 
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print*  ,f  name 

open  ( 1 ,  f  ile=>f n^une ,  f orni»"vinf ormatted" ) 
write (1)  datinl 
close (1) 
else 

pause  '  wrong  integer  in  ifield,  try  again  ' 
goto  1 
end  if 
c 

999  end 

Subroutine  sreinfld(iduin,teinp) 

c************************************************************* 

c  Preset  Fourier  vector  temp  with  reindom  phase  and 

c  unity  euiplitude  coefficients, 

c 

c  Rolcind  Heinrich  Cambridge,  July  24,  1991 

c****  ********************************************************* 

c 

Double  complex  chelp, cim 
Double  Precision  teiip,rauTd,pi 

parameter  (ndim  =  2,itmax  =512,kzmax=128,ivec=itmax*kzmax*2) 
Dimension  temp(ivec) 
c 

pi  =  4 .dO*datan(l.dO) 
cim  =CO.dO,  l.dO) 
c 

do  ii=l,ivec,2 

rand  =  dran2  ( idum) 

chelp  =  2exp(cim*rand*2 .dO*pi) 

temp(ii)  =  dble(chelp) 

temp(ii+l)  =  dimag  (chelp) 

end  do 
c 

return 

end 

Double  Precision  Function  dPAN2(idum) 

c**  *********************************************************** 

c 

c  Rcindom  Function  generator,  see  Numerical  Recipes,  Press  et  al. 

c  page  192.  (Made  double  precision  by  Roland  Heinrich) 

c  Returns  a  uniform  remdom  deviate  between  0.0  and  1.0.  Set 

c  "idum"  to  any  negative  value  to  initialize  or  reinitialize 

c  the  sequence, 

c 

c  Roland  Heinrich  Ccunbridge,  Feb  5,  1991 

C 

Dimension  ir(97) 
double  precision  rm 

Parameter (m=7 14 02 5,  ia  =  1366,  ic=l 50889, rm=l.dO/M) 

Data  iff  /O/ 
c 

if (idum.lt.O.or.iff .eq.O)  then 
iff  =1 

idum  =  mod  (ic- idum,  m) 
do  11  j=l,97 

idum  =  iDod(ia*idu]iH-ic,m) 
ir(j)  =  idum 
11  continue 

idum  «  mod(ia*idu]iH-ic,m) 
iy  =  idum 
endif 

j  =  1  +  (97*iy)/m 
if(j.gt.97.or.j.lt.l)  pause 
iy  =  ir(j) 
dran2  =  iy*rm 
idum  =  mod(ia*idum+ic,m) 
ir(j)=  idum 
c 


return 
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PROGRAM  dstwpO 

C  ********************************************  ...itt  *****************  ******* 

C 

C  stability  progran  for  parallel  three-dixaensional  Blasius  boundary 
c  layer  (spatial  stability) .  Ibis  program  is  based  on  MiXe  Casters 
c  rapid  series  expansion  eigenvalue  calculation  using  Shanks 
c  tremsformation  to  accelerate  the  convergence  process  (qqsum) . 
c 

c  This  program  calls  the  following  subroutines 

c  qnewt  -  multi-dimensional  newton  raphson  root  finder 

c  ludcmp -  l.u.  decomposition  for  qnewt 

c  lubksb -  used  with  ludcnp  to  solve  linear  system  of  eqts. 

c  userfun -  provides  function  and  derivatives  for  qnewt 

c  qfsub  -  rapid  series  eigenvalue  expemsion 

c 

c  The  error-integer  ''ierr''  should  be  zero  for  normal  performance, 

c  It  is  given  the  value  1000  if  the  iteration  iter  in  this  main 

c  program  part  fails  to  converge  within  tolereuice  tol. 

c  For  each  call  to  the  newton  subroutine  qnewt  which  leads 

c  to  non-convergence  within  the  subroutine  in  ntrial  attempts, 

c  ierr  is  incremented  by  1. 

c 

c  The  first  three  sx±)routines  are  adapted  from  "Numerical  Recipes" 

c  William  H.  Press,  Briein  P.  Flannery,  Saul  A.  Teukolsky, 

c  William  T.  Vetterling,  Cambridge  University  Press,  1986. 

c 

c  October  9,  1990  modified  to  double  precision  emd  use  of 

c  qfsub  (Roland  Heinrich) 

c  October  17,  1990  modified  for  "B"  version,  i.e.  outer  loop  is 

c  eliminated  and  dfuserb  is  used  (hch) . 

c  October  18,  1990  this  program  dst_wpB  is  modified  dstabldB 

c  to  allow  mapping  of  Re,  spanwise  wavenumber  br 

c  space.  The  physical  experimental  parameters  are 

c  set  in  DATA  statement.  Make  toleramces  of  convergence 

c  test  value -dependent  (atol,  omega  tol) 

c 

c  November  6,  1990  extrapolate  guess  for  alpha  from  previous 

c  calculated  values  (Hch) . 

c 

c  November  10,  1990  preset  all  eigenvalues  "ca"  to  (1.0, 1.0) 

c  and  calculate  only  up  to  frequncy  (  in  frequency 

c  loop  Iomega)  where  tolerance  criterium  is  missed  for 

c  the  first  time.  Also  "taper  off"  the  last  (in  omega) 

c  value  s.t.  a  smooth  frequency  window  is  used  — >  avoids 

c  ringing.  At  same  time  set  ierr  to  -99  as  flag, 

c  Read  runpareuneter  from  file  HeadwpD.  (Hch) 

c 

C  ROLAND  HEINRICH  CAMBRIDGE,  MAY  29,1990 

C  *********************************************************************** 

c 

parameter  (np=2) 
c 

implicit  real *8  (a-h,o-y) 

double  precision  rel(32,32),  img(32, 32) , eta, delta, dnu 
double  precision  arcenr,arceni,2u:Tad,ascenr,asccni,asrad 
double  precision  pomst,pomdel,pbdel 
double  complex  cim,ca,castrt,castrtl, chip, chlpl, can, caml 
1  ,castrtlml,castrtml 
double  complex  ccalp(9,20,55) 

dimension  rRe(9,20,55) ,rbr(9,20,55) ,romr(9,20,55) 
dimension  ierarr(9,20,55) 
c 

dimension  x(np) 

C 

common  /inp2/  Rehat, Reihat,omsetr,omseti,br 
common  /head04/  eta, delta, epsO,epsl,arcenr,arceni, 

1  arrad,ascenr,asceni,asrad,m,  im,  istep 
common  /coef04/  rel,img 
c 

c  ******  Thf.  following  data  statement  centaino  tiararvt-  frr  * 


c  performed  augl4,  90,  UO  =  20  m/s,  nu  *  l.57e-5  in*ni/6.  VariadDles  starting 

c  sith  P  are  Physical  varlalbes  %fhlch  have  to  be  nondimens ionali zed 

c  for  run  using  this  program, 

c 

DATA  U0,dnu  /20.d0,1.57d-5/ 

DATA  Restart ,  Redel ,  Hostart ,  castrt / 1260 .  dO ,  65 .  dO , 

1  0.d0,(.15d0,.005d0)/ 

c 

c  ******  read  in  run-parameters 
c 

open(5,file='HeadwpD' , fonii=' formatted') 

read(5,950)  Ponst 

read(5,950)  Poindel 

read (5, 950)  Pbdel 

read (5, 960)  ionax 

read (5, 960)  ibmax 

read (5, 960)  itennax 

close (5) 

950  format(10X,D10.4) 

960  format (lOX, 15) 
c 

print*, 'Pomst  =  ',Pomst,'  Poindel  =  ',Pomdel,'  Pbdel  =',Rx3el, 

1  '  ionax  =  ',  icmax, '  ibnax  =  ', ibmax,'  itermax  =  ',itermc»x 


cim 

dum 

PRe 

castrtl 

Ponur 

tolfac 

iRenax 


(0.d0,l.d0) 

dnu/UO/UO 

Restart 

castrt 

Pomst 

l.d-4 

9 

0 


c****  Preset  all  eigenvalues  with  (l.,1.0) 
c 

do  i  =  1, iRenax 
do  ii  =  l,iboax 
do  iii  =  1,  iomax 

ccalp(i,ii,iii)  =  (I.d0,l.d0) 
ierarr(i,ii,iii)  =  -99 
end  do 
end  do 
end  do 


c  ***  %«rite  Header  in  eignap 
c 

open  (10,  file  =  'Head04' ,  form  =' formatted' ) 
read  (10,*)  m,eta,im 
read  (10,*)  delta,epsO,epsl,istep 
read  (10,*)  arcenr,2u:ceni,arrad 
read  (10,*)  ascenr,asceni,asrad 
close (10) 
c 

open  (10, file  =  'fcoef',form  =  'formatted') 
read  (10,15)  ((rel(i,  j)  ,i=l,iin) ,  j=l,im) , 

1  ((img(i,j)  ,i=l,iin),j=l,im) 
close (10) 

15  format (2X,E24.16E3) 
c 

open(13,file  =  'eigmap2' ,form  =  'formatted') 
c 

write  (13,300)  Restart ,  Redel ,  Pomst ,  Poindel ,  Pbstart ,  Pbdel 
300  format (2X, 'restart  =  ',f6.0,2X,'  delre  *=  ',F6.0,/,2X 

2  ,'  PcOTstart  =  ' 

1  ,f8.3,2X,'  Pdelomr  =  '  ,f6.3,/,2X, '  Pbstart  =  ',F6.3, 

3  2X,'  Pbdel  =  ',f8.3,/  ) 
c 

c  ***********  OUTER  loop  490  in  Re  'iren  '***************** 


do  490  iron  -  1, iRenax 


o  o 
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Rehat  =  Pre 
Reihat  =  O.dO 
castrt  =  castrtl 
chlpl  =  castrtl 
c 

c  Progreun  parameter  Re  dependent 

bst  =  Pre  *  Pb  *  dum*UO 
bdel  =  Pre  *  Ptxlel  *  dun*UO 
ojnrdel  =  Pre  *  Pomdel*  dum 
omstrt  =  Pre  *  Pomr  *duin 
c 

c  keep  loop  st^u^ting  value 

br  =  bst 
c 

c  ************ioop  550  in  s()anwise  wavenumber  'ib'  ******** 
c 

do  550  ib  =  1,  ibmax 
ca  =  castrt 
carol  =  castrt 
chip  =  castrt 

atol  =  zabs  (castrt)  *tolfac 

keep  loop  starting  value 
omr  =  omstrt 
c 

c  ************  inner  frequency  loop  'iomega'  ****************** 
c 

do  551  iomeg  =  l,iomax 
omsetr  =  omr 
omseti  =  O.dO 
omtol  =  omsetr*tolfac 
c 

c  ******  start  Newton-Raphson*  with  3-D  pararoters***** 
c 

ierr  =  0 

can  =  2.d0*ca  -  carol 

x(l)  =  dble(can) 

x(2)  =  dimag(ccin) 

c 

call  dqnewt(itermax,x,2,atol,orotol,ierr) 
c 

caual  =  ca 

ca  =  dciq)lx(x(l)  ,x(2) ) 

c 

write (13, 900)  PRe,amr, test4 
write (13, 910)  ca,ierr 
c 

ccalp(iren,ib,  iomeg)  =  ca 
rbr(iren,ib,ioroeg)  =  br 
romr(iren,ib,  iomeg)  =  omr 
RRe ( iren , ib , iomeg)  =  PRe 
ierarr(iren,ib,  iomeg  )=  ierr 
c 

if  (iomeg. eq.l)  then 
castrtml  =  chip 
chip  =  ca 

ceistrt  =  2.d0*chlp  -  castrtml 
if  (ib.eq.l)  then 
castrtlml  =  chlpl 
chlpl  =  ca 

castrbl  =  2.d0*ca  -  castrtlml 
end  if 
end  if 
c 

c  *****  if  ierr  =»  1,  end  omega  loop  and  jump  to  next  "b"  value 
c  since  series  is  starting  to  have  convergence  prc^lems. 

c  remember  the  higher  omega  values  lead  to  stronger  decaying 

c  eigenvalues  alpha  (provided  I  eun  above  steibility  "bemama"  as 

c  has  to  be  the  case  for  this  exit  procedure  to  sbe  senes ible)  . 
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c 

if  (ierr.eq.l)  then 
iomfin  =  iomega 
goto  552 
end  if 
c 


C 


c 

c 


551 


*****  Update  omega 

omr  =  onr  +  omrdel 
continue 


c 


iomfin  =  iomax 


c 

c  ***  Window  off  eigenvalues  smoothly  with  cos(45,90, 135) . 
c 

552  ccalp(iren,ib,iomfin-2)  =  .85355339d0*ccalp(iren,  ib,  iomfin-2) 

ccalp(iren,ib,iomfin-l)  =  .5d0*ccalp(iren,ib,iomfin-l) 
ccalp(iren,ib,  iomfin)  =  .14644660d0*ccalp(ircn,ib,  iomfin) 
c 

c  *****  update  spanwise  wavenumber 
c 

br  =  br  +  bdel 


write(13,*)  '**************  br  =  ' ,br, ' *******' 
550  continue 


c 

c  *****  update  physical  Reynoldsnumber 
c 

PRe  =  PRe  +  Redel 

write(13,*)  •rnmmm  PRe  *  ^PRe, 
490  continue 


c 

open  ( 1 ,  f  ile="eigwp'' ,  form=''unformatted" ) 
write (1)  ccalp 
write (1)  romr 
write (1)  rbr 
%#rite(l)  RRe 
write (1)  ierarr 
close (1) 
c 

900  fonnat(3x,F6.0,3X,F6.3,3X,E12.3) 

910  fonnat(3x,2F9.6,3X,I6) 

close (13) 


SUBROUTINE  d^lEWr(NTRIAL,X,N,TOLX,TOU',  ierr) 

C  *************************************************************************** 

C  Given  an  initial  ^ess  X  for  a  root  in  N  dimension,  take  NTRIAL  Newton- 

C  R2dipson  steps  to  improve  the  root.  Stop  if  the  root  converges  in  either  sum 

C  vzu:izd>le  increments  TOIX  or  summed  function  values  TOIF. 

C 

C  INPlfT  :  A(i,j)  =  del (fi) /del(xj) 

C  B(i)  =  -  fi 

C  Ref:  Numerical  Recipes,  W.H.  Press  et  al.  pg.  273 
c 

c  version  for  eigenvalue  calculation  gaster  series 

c  ierr  is  error  indicator.  Should  be  zero  for  normal  run.  Each  time  the 
c  tolerance  criterium  is  not  met  for  ntrial  attests,  ierr  is  incremented 
c  by  one. 

C 

c  September  24,1990  modified  for  double  precision  (Hch) 

c 

C  Rolemd  Heinrich  Ceunbridge,  May  29,  1990 

Q  *************************************************************************** 

C 


PARAMETER  (NP  =  2) 


oooooooooooooooo  o  o  non 


Jra* 


implicit  real*8  (a-h,o-z) 

DIMENSION  X(NP),  ALPHA(NP,NP) ,  BETA(NP) ,  INDX(NP) , xinit (NP) 
c 

common  /inp2/  Rehat, Reihat,oinsetr,oinseti,br 
C 

c  *****  keep  initial  guess  for  error  message 
c 

do  2  i  =*  l,np 
xinit (i)  =  x(i) 

2  continue 
c 

DO  13  K  =  1,NTRIAL 

CALL  dUSERFU(X, ALPHA,  BETA) 

ERRF  =  O.dO 
DO  11  1=1, N 

ERRF  =  ERRF+ABS(BETA(I)) 

11  CONTINUE 

IF  (ERRF.LE.TOLF)  RETURN 
C 

CALL  dLUDCMP(ALF«A,N,NP,INDX,D) 

CALL  dLUBKSB(ALF«A,N,NP,INDX,BETA) 

ERRX  =  O.dO 
C 

DO  12  I  =  1,N 

ERRX  =  ERRX  +  ABS(BETA(I)) 

X(I)  =  X(I)  +  BETA(I) 

12  CONTINUE 

IF  (ERRX.LE.TOLX)  RETURN 

13  CONTINUE 


******  write  message  if  tolerance  criteria  not  met  *** 


300 


print* , 
print* , 
print* , 
print* , 
print* , 
print* , 


Tolerauicc  criteria  in  qnewt  not' 
met  within  NTRIAL  attempts  for  case' 

Rehat  =  ', rehat,'  omsetr  =  ',omsetr,'  br  =  ',br 
ntrial  =  ',  ntrial 
errf  =  ',err£,'  tolf  =  ',tolf 

errx  =  ',errx,'  tolx  =  ',tolx 

print*,'  The  initial  guess  for  ca  was  ' ,xinit(l) ,xinit(2) 
print*,  '  Ihe  failed  ca  value  =  ',x(l),x(2) 
print*, 'ca  value  reset  to  initial  guess  before  return' 
do  800  i  =  l,np 
x(i)  =  xinit (i) 
continue 

ierr  =  ierr  +  1 


RETURN 

END 


SUBROUTINE  dUSERFU(X, ALPHA, BETA) 

************************************************************************** 
User  supplied  function  2Uid  derivative  for  multiple  dimension2d 
newton'-rahpson  according  to  NUMERICAL  RECIPES  by  W.H.  Press  p.  271. 

this  program  is  used  for  eigenvalue  calculation  Caster  series 


$$$$$ 

$$$$$ 

$$$$$ 

$$$$$ 


September  22,1990  modified  to  double  precision  (Hdi) 

October  9,  1990  modified  to  use  dqfsum  (Hcd^) 

October  16,  1990  modified  to  take  variation  of  Re(2-d)  into 

acount  when  calculation  domega/dx.  Here 


X  is  3d  wavenumber  "a”  and  ALPHA 
contains  derivatives  dootega/da  (Hch) 


Roland  Heinrich  Cambridge,  May  29,  1990 

************************************************************************** 


PARAMETER  (NI^2) 

C 

implicit  real*8  (a-h,o-z) 

double  precision  dxr,dxl,drer,drci. 


I 


~'i0‘ 


1  delxr 

coinplex*16  ca, cal , com, coml , calp, calpl , cRe, cRel 
c 

DIMENSION  X(NP) ,  AIJ>HA(NP,NP)  ,  BETA(NP) 
c 

common  /inp2/  Rehat, Reihat,omsetr,oinseti,br 
c 

delxr  =  .0005d0 
c 

DO  20  J  =  1,NP 
DO  10  I  =  1,NP 
ALPHA(I,J)  =  0. 

10  CONTINUE 

BETA(J)  =  0. 

20  CONTINUE 

c 

c  ******  calculate  corresponding  2-d  parameters 
c  according  to  parallel  flow  Squire's  transform 

c 

ca  =  dcitplx(x(l)  ,x(2) ) 

calp  =  cdsqrt(ca*ca  +  br*br) 

ere  =  Rehat*ca/calp 

c 

cc  print*,'  2d-parameters  ere, cbet, calp  ' 

cc  print*,  'ere  =  ',cre 

cc  print*,  'cbet=  ',cbet 

cc  print*,  'calp=  ',calp 

c 

drer  =  dble(cre) 

drei  =  dimag (ere) 

dxr  =  dble(calp) 

dxi  =  dimag  (calp) 

C 

c  ****  note  IQ  has  to  be  14,18,22, .  for  subroutine 

c  dqfsum  which  uses  full  shanks 

c 

IQ  =  22 

call  dqfsum(dxr,dxi,drer,drei,IQ, l,dbetr,dbeti) 
c 

c  ****  calculate  3-d  omega 
c 

com  =  ca/calp*dcmplx(dbetr,dbeti) 
c 

BETA(l)  =  -(  dble(com)  -  omsetr  ) 

BETA(2)  =  -(  diroag(ccan)  -  omseti  ) 
c 

cal  =  dcrplx(x(l)+delxr,x(2) ) 

calpl  =  odsqrt(cal*cal  +  br*br) 

cRel  =  Rehat*cal/calpl 

c 

dRerl  =  dble(cRel) 

dReil  =  dimag  (cRel) 

dxrl  =  dble  (calpl) 

dxil  =  dimag  (calpl) 


call  dqfsum(dxrl,dxil,dRerl,dReil,IQ, l,dbetrl,dbetil) 
c 

c  ****  calculate  3-d  omega 
c 

coml  =  cal/calpl*dcmplx(dbetrl,dbetil) 
c 

fprr  =  (dble (coml  -  com))  /delxr 
fpri  =  (dble(coml  -  con))  /delxr 
c 

cc  afprr  =  dabs (fprr) 

cc  afpri  =  dabs(^ri} 

cc  if((afprr.lt.l.d-8)  .or.  (afpri.lt.l.d-8))  then 

cc  print*, 'pause  in  duserfun.f,  slope  small  ' 

cc  print*,'  fprr  =  ',fprr,'  fpri=  ',fpri 


-li¬ 


ce  pause 

cc  else  if  ( (afprr.gt.  I.d5)  .or.  (afpri.gt .  I.d5) )  then 

cc  print*, 'pause  in  duserfun.f,  slope  large  ' 

cc  print*,'  fprr  =  ',fprr,'  fpri=  ' ,  fpri 

cc  pause 

cc  end  if 

c 

ALPHA(1,1)  =  fprr 
alpha (1,2)  =  -fpri 
alpha(2,l)  =  -alpha(l,2) 
alpha(2,2)  =  alpha(l,l) 

C 

RETOPN 

END 

c 

c  *************************************************************************** 

C  Program  QFSum 

C  *************************************************************************** 

C 

C  Fast  version  of  series  program  Feb  1990 
c  uses  Quick  summation  but  Full  matrix  shanks 
C 

SUBROUTINE  DQFSUM  (ATPHAR, ALPHAI,REYR,REYI, IQ, INUM,OMEGaR,CMEGaI) 

C 

implicit  real*8  (a-h,o-z) 

Double  Precision  Real (32,32) ,  Imag(32,32),  Cr(32),  Ci(32) 

Double  Precision  Omegar,  Om^ai,  Realsum,  Imagsum 
Double  Precision  Eta,  Delta,  Br,  Bi 
Double  Precision  ARcenr,  ARceni,  ARrad 
Double  Precision  AScenr,  ASceni,  ASrad 
Double  Precision  Alphar,  Alphai,  Reyr,  Reyi 
Double  Precision  Ur,  Ui,  Vr,  Vi,  Zr,  Zi,  Temp 
Double  Precision  Clr,  Cli,  Ctr,  Cti,  Clmag,  CCr,  CCi 
c 

common  /head04/  eta, delta, epsO,epsl,arcenr,arceni, 

1  arrad,ascenr,asceni,asrad,m,  im,  istep 
common  /coef04/  real,imag 
C 

Integer  R 
C 

Tol=l.d-10 

C 

IR=IQ+1 

C 

Do  100  In=l,Inum 

Ur= Alphar*  Alphar-Alphai*Alphai-AScenr 
Ui=2 .0D0*Alphar*Alphai-ASceni 
Vr^Alphar*Reyr-Alpha  i  *Rey  i  -ARcenr 
Vi=Alphar*Reyi+Alphai*Reyr-ARceni 
Ur=Ur/ASrad 
Ui=Ui/ASrad 
Vr=Vr/ARrad 
Vi=Vi/ARrad 
C 

Ten?>=4Jr*Ur+Ui*Ui 

Br=Ur/TeiBp 

Bi=-Ui/Teinp 

C 

Tenp=Vr*Vr+Vi*Vi 
Zir=(Ur*VrfUi*Vi)  /Temp 
Zi=(-Ur*Vi+Ui*Vr)  /Temp 

C 

C  Diagonal  summation  of  series 
C 

CX:r=0.0D0 

cx:i=0.0D0 

c 


c 


Do  20  R=1,IR 


-  5i- 

Te]np=Br*Vr-Bi*Vi 

Bi=Br*Vi+Bi*Vr 

Br=Teinp 

C 

Realsu]n=0 .  ODO 
Inagsum^O . ODO 
C 

Do  10  J=R,1,-1 
I=R-J+1 
C 

RealsunpRealsunH-Real  (I ,  J) 

Iinagsuin=ImagsunH-Iinag  ( I ,  J ) 

C 

Teiiip=Rea  1  sum*  Z  r- Imagsum*  Z  i 
Imag  suin=Real  sum*  Z  i + Imagsum*  Z  r 
Realsum=Temp 
C 

10  Continue 

C 

CCr=CCr+Br*Realsum-Bi*  Imagsum 
CCi=CCi+Br*Imagsum+Bi*Reals\am 
C 

C  Write  (*,*)  R,  ccr,  CCi 

Cr (R) =ccr 
Ci(R)=CCi 
C 

20  Continue 

C 

C  Shanks  Transformation  —  Full  Matrix 
c 

L  =  0 
M  =  IQ 
c 

40  Do  50  J=2,M 

Jmiv=J-l 
Jpl  =J+1 

Clr=Cr  ( Jmn)  +Cr  (Jpl )  -Cr  ( J )  -Cr  ( J) 

Cl  i=ci  (Jmn)  +Ci  (Jpl )  -ci  ( J )  -Ci  ( J) 

Clmag=Clr*Clr+Cli*Cli 
If  (Clmag  .Lt.  Tol)  Then 
Cr (Jmn)=Cr(J) 

Ci (Jmn)=Ci  (J) 

Else 

ctr^cr  (Jmn)  *Cr  (Jpl )  -ci  (Jmn)  *ci  (Jpl)  -Cr  ( J)  *Cr  ( J)  -K:i  ( j)  *ci  ( J) 
Cti=Cr  (Jmn)  *Ci  (Jpl )  +ci  (Jmn)  *Cr  (Jpl )  -2 . 0D0*Cr  ( J)  *Ci  ( J) 
Cr(Jmn)  =  (Ctr*Clr+Cti*Cli)  /Clmag 
Ci  (Jmn)  =  {-Ctr*Cli+Cti*Clr)  /Clmag 
End  If 

50  Continue 
C 

M  =  M-2 

if  (M.Ne.O)  goto  40 
if  (L.eq.l)  goto  80 
L=1 
c 

M  =  IQ/2  -1 
IM=  M+1 
c 

Do  70  J=2,IM 
Cr(J)  =  Cr(2*J-l) 

70  Ci(J)  =  Ci(2*J-l) 

Goto  40 
c 

80  Omegar  =  Cr(l)*Alphar  -  Ci(l)*Alphai 

Oroegai  =  Cr(l)*Alphai  +  Ci(l)*Altrfiar 
100  Continue 
return 
End 
c 


I 


SUBRCTJTINE  dUJBKSB (A, N , NP, INDX, B) 


n  o  o  o  n  o  o  o  o  o  o  o  o  on  o  o  o  n  o  o  o  n  n  n  o  n  o  o  n  n  o 
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Solves  the  set  of  N  linear  equations  A.X  «  B.  Here  A  is  input,  nat  as  the 
but  rather  as  its  "LU"  decomposition,  determined  by  the  routine  ''UJDCMP. 
is  input  as  the  permutation  vector  returned  by  LUDCMP.  B  is  input  as  the  r 
side  vector  B  and  returns  with  the  solution  vector  X.  A,  N,  NP  and  INDX  ar 
modified  by  this  routine  euid  Ccin  be  left  in  place  for  successive  calls  wit 
different  right-hand  sides  B.  This  routine  takes  into  account  the  possibli 
that  B  will  begin  with  many  zero  elements,  so  it  is  efficient  for  use  in 
matrix  inversion 

Ref.:  Numerical  Recipes,  W.H.  Press  et  al.  pg  37 
September  22,  1990  modified  for  double  precision  (Hch) 

Roland  Heinrich  Cambridge,  May  29,  199 

******  til****  ******  It******************************************************* 

implicit  real*8  (a-h,o-z) 

DIMENSION  A (NP,NP) ,  INDX(N),  B(N) 

II  =  0 

DO  12  1=1, N 

LL  =  INDX(I) 

SUM  =  B(LL) 

B(LL)  =  B(I) 

IF  (II.NE.O)  THEN 
DO  11  J=II,I-1 

SUM  =  SUM  -  A(I,J)*B(J) 

11  CONTINUE 

ELSE  IF  (SUM.NE.O.)  THEN 
II  =  I 
ENDIF 

B(I)  =  SUM 

12  CONTINUE 

DO  14  I=N,1,-1 

SUM  =  B(I) 

IF(I.LT.N)  IHEN 
DO  13  J=I+1,N 

SUM  =  SUM-A(I,J)*B(J) 

13  CCWTINUE 
ENDIF 

B(I)  =  SUM/A(I,I) 

14  CONTINUE 

RETORN 

END 

SUBROUTINE  dLUDCMP(A,N,NP,INDX,D) 

************************************************************************** 
Given  N  xN  matrix  A,  with  physical  dimension  NP,  this  routine  replaces 

"LU"  deconposition  of  a  rowwise  permutation  of  itself.  A  and  N  are  irput. 

,  arranged  as  in  equation  (2.3.14)  of  "Numerical  Recepices  "  W.  H.  Press  e 
pg.  35.;  INDX  is  an  output  vector  which  records  the  row  permutation  effect 
the  partial  pivoting;  D  is  output  as  +/-  1  depending  on  whether  the  number 
interchanges  was  even  or  odd,  respectively.  This  routine  is  used  in  combin 
with  LUBKSB  to  solve  linear  equations  or  invert  a  matrix. 

September  22,  1990  modified  for  double  pecision  (Hch) 

Roland  Heinrich  Cambridge,  May  29,  1990 

************************************************************************** 


PARAMETER  (NMAX=100,  TINY=1.0d-20) 
implicit  real*8  (a-h,o-z) 

DIMENSION  A(NP,NP)  ,INDX(N)  ,W(NMAX) 

D  =  1. 

DO  12  I  =  1,N 
AAMAX=0.d0 
DO  11  J  =  1,N 

IF  (ABS{A(I,J))  .GT.AAMAX)  AAMAX=ABS (A (I ,  J) ) 
11  CONTINUE 


c 


12 


IF  (AAMMC.EQ.  0. )  then 

PAUSE  'SINC3ULAR  MATRIX.  ' 

end  if 

VV{I)  =  l./(AAMAX) 

CONTINUE 


13 

14 


15 


16 


17 


C 


C 


DO  19  J=1,N 

DO  14  I»1,J-1 
SUM=A(I,J) 

DO  13  K«1,I-1 
SUM  =  SUM  - 
CONTINUE 


A(I,K)*A(K,J) 


A(I,J)  =  SUM 
CONTINUE 
AAMAX  =  0. 


DO  16  I=J,N 

SUM  =  A(I,J) 

DO  15  K=l,J-l 

SUM  =  SUM-A(I,K)*A(K,J) 

CONTINUE  '  '  ^ 

A(I,J)  =s  SUM 

DDM  =  W(I)*ABS(SUM) 

IF  (DUM.GT. AAMAX)  THEN 
IMAX  «  I 


AAMAX  =  DUM 
ENDIF 
CONTINUE 

IF  (J.NE. IMAX) THEN 
DO  17  K=1,N 
DUM=A(IMAX,K) 
A(IMAX,K)  =  A(J,K) 
A(J,K)  =  DUM 
CONTINUE 
D=  -D 


W(IMAX)  =  W(J) 
ENDIF 


INDX(J)=IMAX 
IF(A(J,J) .BQ.O.)  THEN 
A(J,J)  *  TINY 


PRINT  *, 
ENDIF 


'A(J,J)  SINGUIAR  (SET  TINY)  FOR  J 


18 

19 


C 


IF(J.NE.N)  THEN 
DUM=  l./A(J,J) 

DO  18  I=J+1,N 
A(I,J)=A(I,J) *D0M 
CONTINUE 
ENDIF 
CONTINUE 
RETURN 

END 


1000,  lO.ODO,  32 
0.000000001,  l.E-20,  l.E-5 
460. ODO,  O.ODO,  255. ODO  ' 
0.063d0,  O.ODO,  0.04275D0 
O.IDO,  -O.ODO 
Poast  =  500. do 
Pondel  »  28. do 

PWel  =  20.  do 

ionax  «  55 

iteax  »  10 

itenoaxs  300 


-‘iS-' 


'Pcoe'f 


Complex  Fourier  coefficients, 
for  wave  packet  eigenvalue  calculation 
i,j,  real  part,  imag.  part 

1  1  0.3223253786563900E+000  0. 1613119803369000E-001 

2  1  -0.2418700978159900E-001  -0.4377522855065800E-003 

3  1  0.7942718453705302E-002  -0. 3595384769141701E-002 

4  1  -0.3036302514374300E-002  0.2014757832512300E-002 

5  1  0.1272812020033600E-002  -0. 1098888693377400E-002 

6  1  -0.5578315467573702E-003  0. 582 15594  61765000E-003 

7  1  0. 252393685514  llOlE-003  -C.3075402346439700E-003 

8  1  -0.1167238078778600E-003  0.1623029093025300E-003 

9  1  0.5490734838531397E-004  -0.8574 1878 13792403E-004 

10  1  -0.2618405960674900E-004  0.4535679545369902E-004 

11  1  0.1263278954866101E-004  -0.2402895006525799E-004 

12  1  -0.6158098585729002E-005  0. 1274826263397699E-004 

13  1  0.3030638026757502E-005  -0. 6772615961381200E-005 

14  1  -0.1505099589849100E-005  0. 3602559672799500E-005 

15  1  0. 754131292524 1098E-006  -0. 1918572706927101E-005 

16  1  -0.3811991859947701E-006  0. 1022883452606000E-005 

17  1  0.1943971028595100E-006  -0.5459203293867201E-006 

18  1  -0.1000188092348300E-006  0.2916541177455700E-006 

19  1  0.5192057628278201E-007  -0. 1559655800065200E-006 

20  1  -0.2719245806304099E-007  0.8348379054723402E-007 

21  1  0.1436661190013000E-007  -0.4472842363156799E-007 

22  1  -0.7655186706756501E-008  0.2398678411452701E-007 

23  1  0.4112458817928700E-008  -0. 1287566142593700E-007 

24  1  -0.2226347994849200E-008  0. 6918002615208295E-008 

25  1  0.1213921518861101E-008  -0.3720592500755500E-008 

26  1  -0. 666228183732 1902E-009  0.2002958021662899E-008 

27  1  0.3677884663488800E-009  -0. 1079366596101700E-008 

28  1  -0.2040865254571101E-009  0.5822549664635800E-009 

29  1  0.1137550123542200E-009  -0. 3144229876461701E-009 

30  1  -0.6364654636659097E-010  0. 1699737572913299E-009 

31  1  0.3572324380751702E-010  -0.9198661971021599E-010 

32  1  -0.2010223064397000E-010  0. 49836884  61912500E-010 

1  2  0.3764735534787199E-001  -0. 1096595078706700E-001 

2  2  •'0.2053064294159401E-002  -0.6840560119599102E-002 

3  2  0.7804390043020202E-003  0.6861797883175301E-003 

4  2  -0.3220866783522100E-003  -0.3602366778068200E-003 

5  2  0.1607301092008100E-003  0. 1076236585504400E-003 

6  2  -0.7966860721353399E-004  -0.4196212466922601E-004 

7  2  0.4099830402992700E-004  0. 1427129063813499E-004 

8  2  -0.2122368641721600E-004  -0.4537243057711700E-005 

9  2  0.1107928710553100E-004  0.9904952094075302E-006 

10  2  -0.5798076472274296E-005  0.1295994422889600E-006 

11  2  0. 303780166177  6000E-005  -0.3866861106871501E-006 

12  2  -0.1590507167748001E-005  0.3618830533014299E-006 

13  2  0.8312615591421498E-006  -0.2705580470774300E-006 

14  2  -0. 433243712905  lOOOE-006  0. 1832859766182100E-006 

15  2  0.2249863797487699E-006  -0.1174688293303900E-006 

16  2  -0.1163220559874400E-006  0.7264628720804501E-007 

17  2  0. 598275207153 1598E-007  -0.4381351459414899E-007 

18  2  -0.3058424624669001E-007  0.2593211156920500E-007 

19  2  0.1552484896194500E-007  -0. 1512300151773600E-007 

20  2  -0.7816018943174198E-008  0.8712952137557302E-008 

21  2  0.3897175915312800E-008  -0.4968416345008102E-008 

22  2  -0. 192098137397 0701E-008  0.2807754917455401E-008 

23  2  0.9337800532804902E-009  -0. 1573949193201000E-008 

24  2  -0.4461136460243198E-009  0.8757866098996203E-009 

25  2  0.2084456218742800E-009  -0 . 4839281197987999E-009 

26  2  -0.9453641341972098E-010  0.2656248832000299E-009 

27  2  0.4109439871724497E-010  -0. 1448558217820300E-009 

28  2  -0.1672276032926900E-010  0.7848858368797504E-010 

29  2  0.6043600949551296E-011  -0.4225309685468299E-010 

30  2  -0.1641670523480600E-011  0.2259558859463200E-010 

31  2  0.1104174362348800E-013  -0. 1200011097435300E-010 

32  2  0.4636372665997600E-012  0.6326588992272598E-011 

1  3  -0.1283412706106900E-001  -0. 3782355459407000E-002 

2  3  0.1200836704811100E-003  0. 5877360235899697E-003 

3  3  -0.2647271030582500E-003  -0. 4775901034008699E-003 
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2 

3 

4 

5 

6 

7 

8 
9 


3  0.1119829903473100E-003 
3  -0.5917461021454098E-004 
3  0.3102428308920901E-004 
3  -0.1607303784 1124 OOE-004 
3  0.8396067642024706E-005 
3  -0.4340464784036202E-005 
3  0.2234226940345300E-005 
3  -0.1138868356065400E-005 
3  0.5746144324803002E-006 
3  -0.2861150107946701E-006 
3  0.14 019794 03 68 13 OOE-006 
3  •'0.6731085733235897E-007 
3  0.3146558213984501E-007 
3  -0.1417664385883200E-007 
3  0. 6046207534  637898E-008 
3  -0.2352392503013301E-008 
3  0.7568264348556397E-009 
3  -0.1219838396959700E-009 
3  -0.9300516606858306E-010 
3  0.1375995567487500E-009 
3  -0.1219906259342100E-009 
3  0.9169635883932202E-010 
3  -0.6338354147095203E-010 
3  0.4161035752070500E-010 
3  -0.2636401340294598E-010 
3  0 . 1626985351999000E-010 
3  -0.9835244597911100E-011 
3  0.5845742726689998E-011 

3  -0.3424968208348700E-011 

4  0.4a91649819910500E-002 
4  -0 . 3271847672294800E-003 
4  0 . 6939046579645904E-004 
4  -0.7333208486670597E-004 
4  0. 28799073334084  99E-004 
4  -0.1663749571889600E-004 
4  0.8228060323745002E-005 
4  -0.4225125849188798E-005 
4  0.2113396703862200E-005 
4  -0 . 10424 1826193 17 OOE-005 
4  0 . 503177432165 lOOOE-006 
4  -0 . 2359472262014601E-006 
4  0 . 1063605097328900E-006 
4  -0.4514414797540701E-007 
4  0 . 1729144294415600E-007 
4  -0. 52910031911324  03E-008 
4  0.5702917027505801E-009 
4  0 . 9721010663099603E-009 
4  -0.1234565671914800E-008 
4  0.1056303933211200E-008 
4  -0.7822532066548095E-009 
4  0.5364329536128300E-009 
4  -0 . 3503781986768901E-009 
4  0 . 22 11493349557600E-009 
4  -0.1360155460483600E-009 
4  0.8194346590162397E-010 
4  -0.4852399732024900E-010 
4  0.2830941731035701E-010 
4  -0 . 1629807053205000E-010 
4  0.9269355248642295E-011 
4  -0.5211722112902198E-011 

4  0. 289803857994  6000E-011 

5  -0.2287265611812500E-002 
5  0. 1306006306549500E-003 
5  -0.8714091381989403E-004 
5  0.2536410647735500E-004 
5  -0.2143877645721700E-004 
5  0.8161139703588599E-005 
5  -0 . 4792507297679500E-005 
5  0.2188780854339700E-005 
5  -0.1083624397324500E-005 


0. 3470053343335201E-004 
-0. 3356829256517800E-004 
0. 3799055093622900E-005 
-0. 2979282101023299E-006 
-0. 1470806296310900E-005 
0. 1396764560013300E-005 
-0. 1073991143130100E-005 
0. 7305732196982700E-006 
-0.4695613426974901E-006 
0. 2903555014 199800E-006 
-0. 1748854998595600E-006 
0. 103248297023 1301E-006 
-0. 5999518037924600E-007 
0.3440246132413398E-007 
-0. 1950165895436798E-007 
0. 1094 1536565667 OOE-007 
-0. 6080690617 693597E-008 
0. 3348950894377400E-008 
-0. 1828337037323999E-008 
0.9895023689310294E-009 
-0.5307969064283201E-009 
0. 2821261002594601E-009 
-0. 1484972145249200E-009 
0.7733847590119001E-010 
-0. 3980798329861202E-010 
0. 2021798180262800E-010 
-0. 1010899176867600E-010 
0. 4959730542930503E-011 
-0. 2376124470793201E-011 
0. 2312668832019000E-002 
-0. 4290250362828401E-003 
0. 4648862159228900E-004 
-0. 5076560773886701E-004 
-0.4128159616811900E-005 
-0. 6333055466711800E-006 
-0. 2947733946712102E-005 
0. 1953499804 585600E-005 
-0. 1497659582128100E-005 
0. 9631414741306800E-006 
-0. 6073055942579197E-006 
0. 3674449828849902E-006 
-0. 2177288394023000E-006 
0. 1264929636590800E-006 
-0. 7236617705075298E-007 
0.4083128501974900E-007 
-0. 2275248611738300E-007 
0. 125293126984 6200E-007 
-0. 6820601861079501E-008 
0. 3670184822723902E-008 
-0. 1951415029566300E-008 
0. 1024401896643200E-008 
-0. 53029164 39298099E-009 
0. 2701995294174300E-009 
-0. 1351464634646901E-009 
0. 6608751740300203E-010 
-0. 3139848797628600E-010 
0. 1434489072665900E-010 
-0. 6186053405632403E-011 
0. 242210179864 0899E-011 
-0. 7745524057203900E-012 
0. 1121861054032500E-012 
-0. 1310433610342400E-002 
0. 1825430226745100E-003 
-0.4431783600011800E-004 
0. 4078320898770499E-006 
-0. 2647962219270999E-005 
-0. 4275434093870000E-005 
0. 1974598490051000E-005 
-0. 1893456555990300E-005 
0. 1127553787228000E-005 


-37- 

10  5  0.4949774847773396E-006  -0.7346052939283297E-006 

11  5  -0.2194350372519700E-006  0. 4371982242901100E-006 

12  5  0.8922443583969692E-007  -0.2590394956314400E-006 

13  5  -0.3135414772259498E-007  0. 1491320773539000E-006 

14  5  0.7247473288174399E-008  -0. 8461537959192398E-007 

15  5  0.1572662S55737699E-008  0. 4720060076124400E-007 

16  5  -0.3923178226727900E-008  -0.2595072423616801E-007 

17  5  0.3820721072855799E-008  0. 1405801075549100E-007 

18  5  -0.2994703152126000E-008  -0.7503086152382801E-008 

19  5  0.2120774666991600E-008  0.3941698967224700E-008 

20  5  -0.1414045547321101E-008  -0.2035103641162600E-008 

21  5  0.9051071558907603E-009  0. 1030052487749300E-008 

22  5  -0.5621230148023002E-009  -0. 5091184 140937299E-009 

23  5  0.3409129367693000E-009  0.2442306912264800E-009 

24  5  -0.2027285977979100E-009  -0. 1125578172334100E-009 

25  5  0.1185277570536800E-009  0.4892324045879800E-010 

26  5  -0.6825583154235203E-010  -0. 1929320245230799E-010 

27  5  0.3875965520760998E-010  0.6211164395308904E-011 

28  5  -0.2171888190127800E-010  -0.9076183814932206E-012 

29  5  0.1201252899235600E-010  -0.9129518972403297E-012 

30  5  -0.6557434283072102E-011  0. 1287445085436400E-011 

31  5  0.3531309142670900E-011  -0. 1139585171020000E-011 

32  5  -0.1874333370802100E-011  0. 86095459 10399905E-012 

1  6  0.1155133359134200E-002  0. 7439779583364700E-003 

2  6  -0.9139500616583998E-004  -0. 8856764 1796544 03E-004 

3  6  0.3818675759248399E-004  0. 1806862383091401E-004 

4  6  -0.2513972867745900E-004  -0. 1255855750059700E-005 

5  6  0.8408450412389399E-005  -0. 3267402917117600E-005 

6  6  -0.6177950126584600E-005  0.2013199491557299E-005 

7  6  0.2259152779515699E-005  -0. 2245070618300801E-005 

8  6  -0.1286934548261299E-005  0. 1237397327713600E-005 

9  6  0.5174964599064003E-006  -0.8772461228545597E-006 

10  6  -0.2316054832363101E-006  0. 5051157359048400E-006 

11  6  0.8239690174605098E-007  -0. 3057180038012999E-006 

12  6  -0.2304902579908199E-007  0. 1740267379091200E-006 

13  6  -0.5660225066428401E-009  -0.9868082884168004E-007 

14  6  0.7497948040224899E-008  0.5451589046856498E-007 

15  6  -0.8083861580132602E-008  -0.2967082757265899E-007 

16  6  0.6574810473836099E-008  0. 1583533659754699E-0C7 

17  6  -0.4747379822589401E-008  -0.8293946862636401E-008 

18  6  0. 32015519 12468400E-008  0.4250661156390799E-008 

19  6  -0.2063437865018600E-008  -0.2124663112113002E-008 

20  6  0.1286505457720500E-008  0. 1029524687723400E-008 

21  6  -0.7814862645893996E-009  -0.4788375806974400E-009 

22  6  0.4645714923867700E-009  0.2099045104397800E-009 

23  6  -0.2710475732747899E-009  -0.8347423446908298E-010 

24  6  0.1554856521313000E-009  0. 2714403528336600E-010 

25  6  -0.8779227345101502E-010  -0.4096140527459802E-011 

26  6  0.4881573270609802E-010  -0. 3905728605119201E-011 

27  6  -0. 267296237227 4299E-010  0.5599025149766200E-011 

28  6  0.1440528512447500E-010  -0.4985926168460402E-011 

29  6  -0.7632294102277701E-011  0. 3780572209455300E-011 

30  6  0.3968000407417202E-011  -0. 2636468300620799E-011 

31  6  -0.2018274653306499E-011  0. 1744281043387900E-011 

32  6  0.9996856857943898E-012  -0. 1111991358368500E-011 

1  7  -0.6207895348779901E-003  -0. 42755494 35988097E-003 

2  7  0.5494606011780000E-004  0.4716699913842600E-004 

3  7  -0.2697741729207300E-004  -0. 6423236300179300E-005 

4  7  0.1194484048028200E-004  -0.5254572101875999E-006 

5  7  -0.7283873856067700E-005  0.2545557663324900E-005 

6  7  0.2561361270636600E-005  -0.2133742100340900E-005 

7  7  -0.1669170615059600E-005  0. 1412515075571700E-005 

8  7  0.5511540734914898E-006  -0. 1007123501040001E-005 

9  7  -0.2873398727842899E-006  0.5764466095570198E-006 

10  7  0.8222573200100701E-007  -0. 3603877303248699E-006 

11  7  -0.2191170089815800E-007  0. 2020414768821900E-006 

12  7  -0.6580677780476696E-008  -0. 1158003541945600E-006 

13  7  0.1250482206671700E-007  0.6329359791834602E-007 

14  7  -0. 1206880551052300E-007  -0. 3429713757441299E-007 

15  7  0.9336679873683801E-008  0. 1805568494717100E-007 


16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


7  -0.6573134481158095E-008 
7  0.4352224802772800E'008 
7  “0.2765316198249899E-008 
7  0.1702247676149700E-008 
7  -0 . 1021468798434500E-008 
7  0.5997095597898300E-009 
7  -0.3452798880143599E-009 
7  0.1951993011672900E-009 
7  -0.1084157416619800E-009 
7  0.5914846390453501E-010 
7  -0.3167016648930598E-010 
7  0.1661336866687301E-010 
7  -0.8513327195980301E-011 
7  0.4241604897176100E-011 
7  -0. 2038861484 157700E-011 
7  0.9328320451154600E-012 

7  -0.3956927559049800E-012 

8  0.3472458338365001E-003 
8  -0.3471606760285799E-004 
8  0.1636026354390200E-004 
8  -0.8201192031265197E-005 
8  0.3631588015196000E-005 
8  -0.2006310069 19 1400E-005 
8  0.6803346082051601E-006 
8  -0.3849222878216101E-006 
8  0.8976159193707606E-007 
8  -0.3217852295733800E-007 
8  -0.1114683900738100E-007 
8  0.1587726572438400E-007 
8  -0.1592112042203600E-007 
8  0.1206583011281700E-007 
8  -0.8507067938978700E-008 
8  0.5597772911869496E-008 
8  -0.3542820481428101E-008 
8  0. 2168131674 196199E-008 
8  -0.1292546403242100E-008 
8  0.7527677925445200E-009 
8  -0.4292299016217500E-009 
8  0.2398201359721200E-009 
8  -0.1312963765487400E-009 
8  0.7036791982439903E-010 
8  -0.3684356639222002E-010 
8  0.1877829662599900E-010 
8  -0.9260928829357701E-011 
8  0.4374123459272001E-011 
8  -0.1941330559848700E-011 
8  0.7774267883099699E-012 
8  -0.2507459300957201E-012 

8  0.3249711358842800E-013 

9  -0.2000474341912200E-003 
9  0.2205145938205500E-004 
9  -0.1035634977597500E-004 
9  0.4981268375558998E-005 
9  -0.2382498905717501E-005 
9  0.1006377033263600E-OOS 
9  -0.4820353751711101E-006 
9  0.13053363 54 562300E-006 
9  -0.5246430845318200E-007 
9  -0.1253537362799800E-007 
9  0.1658314019437100E-007 
9  -0.1960625617414300E-007 
9  0.1456331055749200E-007 
9  -0.1056949905375900E-007 
9  0.693486246206C201E-008 
9  -0.4414659304785600E-008 
9  0.2696528556000999E-008 
9  -0 . 1604856136872000E-008 
9  0.9302206227523900E-009 
9  -0.5269673031271800E-009 
9  0.2917187880147099E-009 


-0. 9305259673908497E-008 
0.4654383989333200E-008 
-0. 22494 30197 665000E-008 
0. 1036899788253700E-008 
-0 . 44  58946822882799E-009 
0. 1699291957146800E-009 
-0.4889092602988798E-010 
0. 1015267514156700E-011 
0. 1419989473228100E-010 
-0. 1608915084494500E-010 
0. 1341851716546398E-010 
-0. 9860862126842601E-011 
0. 6746944582885996E-011 
-0.4403598145852301E-011 
0. 2776854050520701E-011 
-0. 1704399859095400E-011 
0. 1022998743008300E-011 
0 . 2501737326383600E-003 
-0. 2471319749020000E-004 
0. 2916585117418400E-005 
0.1808303068173700E-005 
-0. 1799942310754001E-005 
0. 1660673888181900E-005 
-0. 1052026163961300E-005 
0. 6777846692784799E-006 
-0. 4142605121160202E-006 
0. 2360507380672100E-006 
-0. 1363840311796600E-006 
0.7404410240496904E-007 
-0.4020103361312998E-007 
0. 2090616568750600E-007 
-0. 1065861177096400E-007 
0. 5213318221564098E-008 
-0. 2442081870057700E-008 
0. 1067627652950900E-008 
-0 . 4 174 0644 2743 1300E-009 
0. 12717724 09605599E-009 
-0. 9771023400106898E-011 
-0. 2930330283978900E-010 
0. 3571273138325298E-010 
-0. 3051425084832401E-010 
0. 2270308940843800E-010 
-0. 1564800025666300E-010 
0. 1025860126013800E-0] 0 
-0. 6484702098213902E-':il 
0.3983525748846300F-011 
-0. 2389480106834699E-011 
0. 1403785182714500E-011 
-0.8091772152503397E-012 
-0. 1484198874095500E-003 
0. 1340333346888700E-004 
-0. 6799712082283800E-006 
-0. 1406387696079000E-005 
0. 1628367613193400E-005 
-0. 1084664631889600E-005 
0.7932791277198704E-006 
-0. 4573246599193199E-006 
0. 2819498376993601E-006 
-0. 1585158457828600E-006 
0. 8800819273346900E-007 
-0.4754302 196374701E-007 
0. 2465320037003900E-007 
-0. 1248735159720100E-007 
0. 6003608721272300E-008 
-0. 2744624527650799E-008 
0. 1142032246548300E-008 
-0. 4032469080872400E-009 
0.8577775295615098E-010 
0. 3235147566504899E-010 
-0. 6345059200274500E-010 


22 

9 

23 

9 

24 

9 

25 

9 
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27 
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28 

9 

29 

9 

30 

9 

31 

9 

32 

9 

1 

10 

2 

10 

3 

10 

4 

10 

5 

10 

6 

10 

7 

10 

8 

10 

9 

10 

10 

10 

11 

10 

12 

10 

13 

10 

14 

10 

15 

10 
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10 
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10 

1 

11 

2 

11 

3 

11 

4 

11 

5 

11 

6 

11 

7 

11 
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11 
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11 

18 

11 

19 

11 

20 

11 

21 

11 

22 

11 

23 

11 

24 

11 

25 

11 

26 

11 

27 

11 

-0. 1577604158420300E-009 
0 . 8317818656067299E-010 
-0.4260441652248502E-010 
0.2106681669389299E-010 
-0 . 9946135928751207E-011 
0 . 4390772034 151800E-011 
-0. 17308656678067 OOE-011 
0 . 5311599684 39 3300E-012 
~0.3949113917280600E-013 
-0 . 12674061832 58100E-012 
0 . 155291822153 1900E-012 
0 . 117852869152 5600E-003 
-0 . 1409151263942500E-004 
0 . 6588044470845498E-005 
-0. 3104950337729000E-005 
0. 1429191115676100E-005 
-0.6154281209092000E-006 
0.2246797379257300E-006 
-0.7614639940811702E-007 
-0.3496972711403098E-008 
0. 15054704 14997500E-007 
-0 . 2234129503619900E-007 
0.1661547699427500E-007 
-0 . 1273098693 13  3099E-007 
0 . 8326819234 127900E-008 
-0 . 5396890490061399E-008 
0. 3294507555651200E-008 
-0. 1969056251383000E-008 
0. 1138652061527500E-008 
-0 . 6431644661297796E-009 
0 . 3536657355862600E-009 
“0 . 1894326639995100E-009 
0 . 9846345266906198E-010 
-0.4941250186796300E-010 
0 . 2370184257194400E-010 
-0.1066802722549100E-010 
0.4329387976592702E-011 
-0. 1417504893845600E-011 
0. 192913989904 1499E-012 
0.2425030391339101E-012 
-0.3360367539784202E-012 
0 . 3002503844 509698E-012 
-0.2295990275073099E-012 
-0 . 7064790406730001E-004 
0.9063448487722797E-005 
-0 . 4191812422504901E-005 
0. 1952019601958499E-005 
-0. 8586345074 945701E-006 
0.3528353147430602E-006 
-0. 1168280192 587200E-006 
0. 2024697209890300E-007 
0. 1179747677326800E-007 
-0.2215465322308300E-007 
0. 1832565210690999E-007 
-0. 1472674604485700E-007 
0.9747264329007496E-008 
-0.6484711878584903E-008 
0 . 3963841255227900E-008 
-0. 2395671661048500E-008 
0 . 13842975654654 OOE-008 
-0. 7832773318838803E-009 
0 . 4288659705142800E-009 
-0. 2283471467468000E-009 
0. 1173495467465800E-009 
-0. 5788050247423899E-010 
0. 2699034919 173600E-010 
-0. 1157052578748500E-010 
0.4249752025980999E-011 
-0. 10210509738054  OOE-011 
-0. 23471326336504  OOE-012 


33- 

0. 6051427414721601E-010 
-0.4737206460436700E-010 
0. 3363663167443901E-010 
-0. 2248764369161500E-010 
0. 1441058210260900E-010 
-0. 8939269463548204E-011 
0. 53992843 18252703E-011 
-0. 3186597036589100E-011 
0. 1841580490533099E-011 
-0. 1043242098611200E-011 
0. 5794199213309001E-012 
0. 8912580960895898E-004 
-0. 7268514764 29 1501E-005 
-0. 8478147606183498E-007 
0. 1268378468921600E-005 
-0. 1119605371968600E-005 
0. 8484583418066904E-006 
-0. 5100789621792501E-006 
0. 3330213758 090399E-006 
-0. 1818749950643900E-006 
0. 1065717682990901E-006 
-0. 5610659670196599E-007 
0. 2977644264490199E-007 
-0. 1488175804809100E-007 
0. 7142134439419600E-008 
-0. 3206127585642100E-008 
0. 1289190865172400E-008 
-0. 4172084344666898E-009 
0. 5202129352288600E-010 
0. 7491031406292603E-010 
-0. 1006242103529500E-009 
0.8851087918 149105E-010 
-0. 6681574737932299E-010 
0. 4641212761336001E-010 
-0. 3054238112420999E-010 
0. 1931929963228099E-010 
-0. 1184318615399200E-010 
0. 7070248900231702E>011 
-0. 4122193905664600E-011 
0. 2350656778601898E-011 
-0. 1311579534216100E-011 
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26  18  -0.9205034737919100E-012  -0.2871745538850900E-012 

27  18  0.5636297635211500E-012  0. 6825832623734102E-013 

28  18  -0.3335330330065500E-012  0. 1980334909186699E-013 

29  18  0. 1910574563 140900E-012  -0.4611789127125501E-013 

30  18  -0.1058666975156500E-012  0. 4631986458346200E-013 

31  18  0.5657082332959098E-013  -0.3756590443446500E-013 

32  18  -0.2895937206400300E-013  0.2741861235668400E-013 

1  19  -0.1665242734816300E-005  -0. 1232271756634900E-005 

2  19  0.2927462219304301E-006  -0. 3818398042199100E-008 

3  19  -0.1202092079211100E-006  0.8056045430748796E-007 

4  19  0.4300697398207395E-007  -0.6893008475117300E-007 

5  19  -0.9695581226765203E-008  0. 4608617487633599E-007 

6  19  -0.2847598157274000E-008  -0. 2740828008995800E-007 

7  19  0.6146634756731797E-008  0. 14 94 98387003 1800E-007 

8  19  -0.5821526816163198E-008  -0. 7491088638289503E-008 

9  19  0.4427766153725100E-008  0. 3372078394292500E-008 

10  19  -0.3008659765768602E-008  -0. 1273625094278000E-008 

11  19  0.1893861512058899E-008  0.2965322987158499E-009 

12  19  -0.1121623238731200E-008  0.9194463940209697E-010 

13  19  0.6291392962154399E-009  -0.2067251092530099E-009 

14  19  -0.3321348196472200E-009  0.2015959621459698E-009 

15  19  0.1646415087597200E-009  -0. 1599562843290800E-009 

16  19  -0.7357072590030702E-010  0. 1129767945862200E-009 

17  19  0,2828085161887800E-010  -0.7455692313529400E-010 

18  19  -0.6830428583848303E-011  0. 4651383792020300E-010 

19  19  -0.168948589927 14 OOE-011  -0. 2771964255354800E-010 

20  19  0.4372076051889501E-011  0. 1583437027330400E-010 

21  19  -0.4396569913689397E-011  -0.8643414109443802E-011 

22  19  0.3567278850785699E-011  0.4497325688940199E-011 

23  19  -0.2583279944123901E-011  -0.2194758264018000E-011 

24  19  0.1748493060407800E-011  0.9809205314334504E-012 

25  19  -0.1123648808547300E-011  -0. 3732985600322898E-012 

26  19  0.6932483558562401E-012  0.9298632149641307E-013 

27  19  -0.4123146865164101E-012  0.2 144 8807 15366502E-013 

28  19  0.2369411633042201E-012  -0.5665280256635798E-013 

29  19  -0.1314502570378200E-012  0.5797878859957100E-013 

30  19  0.7016944263051704E-013  -0.4738926400816102E-013 

31  19  -0.3578055549329600E-013  0. 3475439201835701E-013 

32  19  0.1718956296476798E-013  -0.2381136161454200E-013 

1  20  0.1069000632014601E-005  0.7812175226717997E-006 

2  20  -0.1917887715308100E-006  0. 1115814374230700E-007 

3  20  0. 7704 101534 500296E-007  -0. 5770498034962700E-007 

4  20  -0.2606804549998300E-007  0.4778972950703099E-007 

5  20  0. 450901005422 llOOE-008  -0. 3138794468782201E-007 

6  20  0. 323644666622 1599E-008  0. 1834421681223800E-007 

7  20  -0. 494512 6086397496E-008  -0.9788228005902502E-008 

8  20  0.4377858076054500E-008  0.4747720883102602E-008 

9  20  -0.3222177857864500E-008  -0.2020841494143600E-008 

10  20  0.2137416243996900E-008  0.6679954922539400E-009 

11  20  -0.1317847275750199E-008  -0.6987165707128497E-010 

12  20  0.7634629040076401E-009  -0. 1487068940209100E-009 

13  20  -0.4168951295291400E-009  0. 1914650382683900E-009 

14  20  0.2131870652233300E-009  -0. 1676847133369500E-009 

15  20  -0. 999954771763 1801E-010  0. 1249230024980000E-009 

16  20  0.4122465563360898E-010  -0. 8549625590825702E-010 

17  20  -0.1259308976542299E-010  0. 5463899124813500E-010 

18  20  0.3672404448682602E-012  -0. 3328030559468499E-010 

19  20  0.4010790831399102E-011  0. 1930373916270200E-010 
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-0.4685853267905001E-011 
0.4032102776663200E>011 
-0.3010285163501699E-011 
0 . 2080630372852700E-011 
-0. 1355666121895200E-011 
0. 84549571120404  04E-012 
-0. 5067386397 647502E-012 
0 . 2929251079978501E-012 
-0.1631322301635201E-012 
0 . 8726899140398105E-013 
-0.4448702419900703E-013 
0.2129047467296599E-013 
-0 . 9252012133850698E-014 
-0.6885908874210100E-006 
0. 1257185004988100E-006 
-0 . 4930735642005901E-007 
0 . 1563355134  237600E--007 
-0.1695700246884300E-008 
-0.3049519747833101E-008 
0.3854779162537598E-008 
-0.3241001023113200E-008 
0.2317770197990400E-008 
-0. 1504570135146801E-008 
0.9080115548343301E-009 
-■0.5143304671051400E-009 
0.2727303938243600E-009 
-0 . 1339577893055799E-009 
0 . 5889875392961502E-010 
-0.2112155242373101E-010 
0.3859537689443497E-011 
0 . 2926236293207500E-011 
-0 . 4676475786474800E-011 
0 . 4391511893714299E-011 
-0 . 3424927008666097E-011 
0 . 2432752567 102501E-011 
-0.1614617836133400E-011 
0 . 1021465030615100E-011 
-0 . 6186817551323596E-012 
0. 3607594908 569498E-012 
-0.2022439561742001E-012 
0 . 1087582104995600E-012 
-0 , 5564386774 156500E-013 
0 . 2667556795029800E-013 
-0 . 1158980273948200E-013 
0 . 41674982 38730698E-014 
0.4449006780760100E-006 
-0.8244310834015797E-007 
0 . 3150394789486200E-007 
-0.9249449206549800E-008 
0.2510442131953300E-009 
0 . 2627291717516300E-008 
-0 . 2937244891754200E-008 
0 . 2368595986013101E-008 
-0. 1651529579760300E-008 
0. 1049622500026699E-008 
-0 . 6203756397482600E-009 
0 . 3426315897670000E-009 
-0 . 1759919293631300E-009 
0.8219799146891299E-010 
-0 . 3305871201786398E-010 
0.9291577056369499E-011 
0 . 8257345003072306E-012 
-0 . 4195587887528801E-011 
0. 4553371135002100E-011 
-0. 3784683504 093300E-011 
0 . 2785616138797899E-011 
-0 . 1894832814958200E-011 
0. 1220405937765700E-011 
-0 . 7495323806061496E-012 
0.4421284031160702E-012 


-0. 1070943420750101E-010 
0. 5646071717796398E-011 
-0. 2798655404 64 5499E-011 
0. 1273324002661100E-011 
-0. 4994330144759999E-012 
0. 1361879359285100E-012 
0. 1580355787097900E-013 
-0. 6526293990904300E-013 
0. 6996603952669600E-013 
“0. 5827106917648304E-013 
0.4317038263326700E-013 
-0 . 2978060099955200E-013 
0. 1951867148268999E-013 
-0. 4964289246345301E-006 
-0. 1286580175730001E-007 
0. 4102103190462002E-007 
-0.3301054007920400E-007 
0. 2131003462579900E-007 
-0 . 1223212997558700E-007 
0.6373102046097800E-008 
-0. 2978435942324599E-008 
0. 1180663011801600E-008 
-0. 3172173912435698E-009 
-0.4411823337213702E-010 
0. 1581772490766300E-009 
-0. 1648709224699700E-009 
0. 1331617593969300E-009 
-0.9549426527311105E-010 
0 . 6309621575217903E-010 
-0 . 3943051440913301E-010 
0 . 23  33907719864699E-010 
-0. 1320400386251000E-010 
0. 7078963717294100E-011 
-0 . 3580538426514701E-011 
0 . 1666729905133699E-011 
-0. 6803554030918002B-012 
0.2061866844795000E-012 
“0. 1695028123936700E-014 
-0. 6970146709 155498E-013 
0. 8128490014327303E-013 
-0. 6983800202638895E-013 
0. 5262553698479202E-013 
-0. 3668024678481602E-013 
0. 2421438836117500E-013 
-0. 1531865353153201E-013 
0. 3161164556786401E-006 
0. 1201823707219800E-007 
-0. 2898052642308400E-007 
0 . 2272755494914200E-007 
-0. 1442504249382600E-007 
0. 8126295192312201E-008 
-0.4124850239151101E-008 
0. 1846066854938000E-008 
-0. 6667474350408000E-009 
0. 1205866240194800E-009 
0.9222711483403103E-010 
-0. 1461854387541100E-009 
0. 1345544092812300E-009 
>0. 1030890442455000E-009 
0. 7121760819961402E-010 
-0 . 4591290195254899E-010 
0. 2789778130729599E-010 
-0. 1613939364097900E-010 
0.8837497574021302E-011 
-0.4574710835841802E-011 
0. 2187910226256200E-011 
-0.9336701333947897E-012 
0. 3136194103893502E-012 
-0. 3716535271810298E-013 
-0. 6682146697946401E-013 
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-0. 2501929598827001E-012 
0. 1356701209709301E-012 
-0. 6091159902510999E-013 
0. 3375463903071197E-013 
-0. 1476593404026300E-013 
0. 5360124016596496E-014 
-0. 1072778027367600E-014 
-0. 2882322576169799E-006 
0. 5407931524814600E-007 
-0. 2008781940787699E-007 
0. 5376358469533200E-008 
0. 4232876826826298E-009 
-0. 2147401589880400E-008 
0. 2199533888358500E-008 
-0. 1712784247764900E-008 
0. 1166811869346400E-008 
-0.7265587576021904E-009 
0. 4200602365944699E-009 
-0. 2259012699 123999E-009 
0. 1116505360387900E-009 
-0.4904103165226402E-010 
0. 1720282249984000E-010 
-0. 2633370735014001E-011 
-0. 3017729295617900E-011 
0. 4397768607611299E-011 
-0. 4025244113720100E-011 
0. 3109438271101101E-011 
-0.2184627045237500E-011 
0. 1438747341750400E-011 
-0.8995635182931798E-012 
0. 5382808998880002E-012 
-0. 3084131799772900E-012 
0. 1691498503764199E-012 
-0. 8811950735071500E-013 
0.4304113218178002E-013 
-0. 1909483991280699E-013 
0. 7108247977315401E-014 
-0. 1579501406 139800E-014 
-0. 6493446159088901E-015 
0. 1871882204795800E-006 
-0. 3547956950455998E-007 
0. 12777752 18966300E-007 
-0. 3051385144559001E-008 
-0.6790498718878302E-009 
0. 1695131701673300E-008 
-0. 16246204 38178700E-008 
0. 1227405399539800E-008 
-0. 8181731403489298E-009 
0.4991507762675001E-009 
-0. 2820302325012799E-009 
0. 1471712474110400E-009 
-0. 6955509779249500E-010 
0. 2809083000932000E-010 
-0. 7842198244956100E-011 
-0. 8770521743757500E-012 
0. 3770331269414799E-011 
-0. 4057162591997198E-011 
0. 3361040178812000E-011 
-0. 2463686589807001E-011 
0. 1669034810532200E-011 
-0 . 1067 1083 1577  34 OOE-011 
0. 6497773239605400E-012 
-0.3781721030077200E-012 
0. 2103391923526399E-012 
-0 . 1111392607718700E-012 
0. 5513771473360398E-013 
-0.2497261668017301E-013 
0.9646056096078194E-014 
-0. 2445232076309100E-014 
-0. 5469981013481796E-015 


0. 9038939301903000E-013 
-0. 8145690881793801E-013 
0. 6288447139315300E-01 3 
-0. 4449877424005100E-0i3 
0. 2967875714610600E-013 
-0. 1891320386166299E-013 
0. 1160716437380200E-013 
-0. 2016701898810399E-006 
-0. 10205860156987 OOE-007 
0. 2036865609511600E-007 
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-0. 1645409503092700E-010 
-0. 1049611786374500E-Oi.  . 
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-0. 62358884 0612 5000E-013 
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0. 1497702481511400E-014 
-0. 1218332101871100E-006 
0.2327819004221999E-007 
-0. 8104988680201998E-008 
0. 1674391625350300E-008 
0.7193800555249696E-009 
-0. 1305358821035400E-008 
0. 1186600373515000E-008 
-0 . 8727025768706696E-009 
0 . 5697650129476999E-009 
-0.3404752590974200E-009 
0. 1876941241318900E-009 
-0.9468035383486397E-010 
0.4233081593585402E-010 
-0. 1519682123196600E-010 
0. 2516799052151900E-011 
0.2478780418771500E-011 
-0 . 3765919433934600E-011 
0.3479556920371602E-011 
-0.2700893328433401E-011 
0. 1898302695591002E-011 
-0. 1247420027515600E-011 
0. 7759768847447100E-012 
-0.4602146029715502E-012 
0.2603005430559100E-012 
-0. 1399135934063500E-012 
0. 7073162178574197E-013 
-0.3284121564105200E-013 
0. 1325090719388100E-013 
-0. 3816422563851598E-014 
-0.2212519541399700E-015 
0. 1610732463817100E-014 
-0. 1805917653984699E-014 
0.7945334346004498E-007 
-0. 1527217996510899E-007 
0.5124197954842201E-008 
-0.8730123290945400E-009 
-0. 6587686507408799E-009 
0.9866780725786800E>009 
-0 . 85859014  34266006E-009 
0. 6162106380713800E-009 
-0. 3942460524708501E-009 
0. 2306178303879101E-009 
-0. 1237790286934300E-009 
0. 600364 5497409897E-010 
-0. 2501605346483200E-010 
0.7457392342535797E-011 
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